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ABSTRACT 

We present photometric and spectroscopic observations of galaxies hosting Type la supernovae 
(SNe la) observed by the Nearby Supernova Factory (SNfactory). Combining GALEX UV data with 
optical and near infrared photometry, we employ stellar population synthesis techniques to measure 
SN la host galaxy stellar masses, star-formation rates (SFRs), and reddening due to dust. We reinforce 
the key role of GALEX UV data in deriving accurate estimates of galaxy SFRs and dust extinction. 
Optical spectra of SN la host galaxies are fitted simultaneously for their stellar continua and emission 
lines fluxes, from which we derive high precision redshifts, gas-phase metallicities, and Ha-based 
SFRs. With these data we show that SN la host galaxies present tight agreement with the fiducial 
galaxy mass-metallicity relation from SDSS for stellar masses log{M^,/MQ) > 8.5 where the relation 
is well-defined. The star-formation activity of SN la host galaxies is consistent with a sample of 
comparable SDSS field galaxies, though this comparison is limited by systematic uncertainties in SFR 
measurements. Our analysis indicates that SN la host galaxies are, on average, typical representatives 
of normal field galaxies. 
Subject headings: supernovae: general 



1. INTRODUCTION 

Type la supernovae (SNe la) are some of the most 
dramatic explosive events in the Universe, yet the ex- 
act nature of their stellar progenitors remains a mys- 
tery. SNe la are very bright, and exhibit low intrin- 
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sic lu minosity dispersion (^ 
dis 



0.35 mag; Branch & Miller 
19931), making them suitable for measuring extragaiac- 
istances. The discovery of an empirical relation- 
ship betwe en the luminosity and light curve decline rate 



of SN e la (Pskovskii 1977 IPhillips 



19991, aided m large part by theaiscovery ot e x treme 



1993 Phillips et al. 



events such as SN 1991T (Phillips ct al. 1992; Filip 
penko et aL][T992aJ and SN 1 991 bg (Filippcnko ct aE 
Tn92b; Lcibundgut et ah 1993), resulted in a significant 



improvement in the standardized luminosity of SNe la 
and enabled their successful use as cosmological probes. 
Discovery of a second empirical relationshi p between the 



observed color and luminosity of SN e la ( Hamuy et al 

further decreased 



1996 Riess et al. 1996 



1998) 



Tripp 

the observed dispersion to the current level of about 
~ 0.15 mag. This uniformity of SN la luminosities has 
enabled their successful use as standardizable candles to 



j.998 
et al 



cons train t he energy content of the Universe ( Riess et al. 
Pcrl mutter et al.|[T999l |Tonry et al.||200i)| jKnop 
Wood-Vasey et al.||2007| |Riess et ar]|2007| 



.2003 



Kowalski et al.||2008l IHicken et al.||2009l IKessler et al. 



2009||Amanullali et al.|2010||Sullivan et al.|2011[|Suzuk 



et al.|2012[ ). 

Despite their utility as cosmological distance indica- 
tors, the stellar progenitors of SNe la remain as yet un- 
determined. With future SN la surveys pl anned to find 
SNe la at high redshifts (e.g. WFIRST - [Green et al 
2012 ), concerns remain that the younger stellar ages and 
lower metallicities of high redshift environments could 
bias cosmological measurements if the corrected bright- 
nesses of SNe la vary with these parameters. A promising 
source for clues to the origin of SN la brightness diversity, 
and its possible dependence on progenitor properties, is 
the study of SN la environments. 

Early studies of SN la host galaxies found qualita- 



tive evidence for a correlation between the observed peak 
magnitude, light curve decline rate, and expansion veloc- 
ity of an SN la wi th the morphological type of its host 



et al. 2010 Lara-Lopez et al 
ical check on this technique s 



2010). However, a crit- 
lould be a confirmation 



galaxy (JFil ippcnko 1989 Branch fc van den Bergh|[l993 
Mamuy et al. 1996), such that brighter slower decimmg 



biNe la preferentially occur in later type (spiral and irreg- 
ular) galaxies, while fainter slower declining SNe la pref- 
erentially occur in earlier type (E/SO) galaxies. Recent 
studies have confirmed these trends and found analogous 
corr elations wit h respect to host galaxy m ass or mctallic- 
ity dHamuy eFal. .2000. Gallaghe r et al.|2005, 2008; NeiU 



et al.||200y[ [Howell et al.||2009] ). 'I'hese treh'Ss'VrOTaHe' 
compelling evidence that some property of the progeni- 
tor (most likely age or metallicity) that correlates with 
host galaxy properties may be driving the observed di- 
versity of SN la light curve decline rates. 

Recently, troubling evidence for progenitor-driven 
trends of corrected SN la luminosities arose from stud- 
ies of SN la Hubble residuals and the properties of their 

120101 




a cor- 

relation can easily bias the measurement of cosmolog- 
ical parameters, particularly the dark energy equation 
of state parameter and its evolution with redshift. This 
subject will be the focus of the second paper in this series, 
but here it motivates a critical question: if a particular 
property of host galaxies should be used to correct SN la 
luminosities, what is the best technique for measuring 
both the value and uncertainty of that property? 

The most accessible physical properties to measure for 
galaxies are those which can be inferred from multi-band 
photometry, particularly stellar mass and star-formation 
rate (SFR). To derive those galaxy properties from pho- 
tometry requires invocation of the art of stellar popu- 
lation synthe sis (SPS). This is a rich field of study (see 
Conroy[2013 for a recent review), to which we will later 



devote an entire section, and is being used with increas 
ing frequency in the study of SN la host galaxies. SPS 
involves the comparison of observed galaxy photometry 
(across multiple filters) to fiuxes predicted for a model 
spectral energy distribution (SED) generated for a cho- 
sen galaxy star- formation history (SFH). While much 
effort in the field of SPS is devoted to deriving model 
SEDs that best match observed galaxy photometry, we 
are here more interested in how much variation of the 
SFH is allowed for a galaxy given its observational mea- 
surement errors and a known level of uncertainties in the 
models. We will also inspect the effects of broader wave- 
length coverage in p hotometry, and will echo the findings 
of previous authors ( Pa povich et al. ' 2001 Gil de Paz fc 



Madore||2002 
showed that 



u 



NeiU et al. 2009 ; Gupta et al. 
V data is critical tor accurate 



i>01l| who 



y assessing 



a galaxy's recent star-formation activity. 

While the photometrically derived properties of stel- 
lar mass and SFR are the most easily accessible galaxy 
properties, metallicity is a property of critical interest for 
SNe la. Galaxy metallicities are most directly measured 
from spectroscopy, but as this is observationally expen- 
sive a more desirable technique would be to infer metal- 
licity from photometrically derived properties. This can 
be achieved by invoking galaxy scaling r elations such as 
the galaxy mass-metallicity relation (e.g. Tremont i et al. 
[2004|) or the mass-metallicity-SFR relation (,Mannucci 



that SN la host galaxies exhibit metallicities and SFRs 
consistent with expectations given their stellar masses, 
which we perform in this work. 

This paper is the first in a series of papers which will 
study the host galaxies of SNe la discovered or followed 
by the Nea rby Supernova Factory (SNfactory, [Aldering 



et al. 2002 1 . Here we present the observational data for 



the SNfactory host galaxy sample which will form the ba- 
sis of the scientific analyses in this and future papers. In 
Section [2] we present the photometric and spectroscopic 
observational data for our SN la host galaxy sample. 
We derive estimates of stellar mass, star-formation rate 
(SFR), and dust extinction for each host galaxy using a 
stellar population synthesis method described in detail 
in Section [3j We then exploit the derived host galaxy 
physical parameters (mass, metallicity, and SFR) to ex- 
amine the question of whether SN la host galaxies are 
typical representatives of the normal galaxy population 
in the local universe, considering SN la host metallicities 
in Section m and host SFRs in Section [H We discuss the 
implications of our analyses in Section |6] and summarize 
our conclusions in Section [T) 

Throughout this paper we employ a standard ACDM 
cosmology with Q\ — 0.7, flM — 0.3, Hq — 
70 km s~^ Mpc""'^. Stellar masses and sta r-formation 



rates in this paper are computed using a Chabrier ( 2003 1 
initial mass function (IMF), and our prescriptions tor 
converting literature values to this choice of IMF are de- 
scribed where appropriate. 

2. HOST GALAXY DATA 

In this Section we present the observational data for 
the full sample of SNfactory SN la host galaxi es. We 
first bri efly describe the SN la sample in Section |2.1[ In 
Section [Z!2| we describe the galaxy photometric data that 
was collected from numerous public sources, as well as 
the reduction of targeted observations using the SNfac- 
tory SNIFS instrument. Derivation of galaxy stellar 
masses and specific star formation rates utilizes a com- 
plex custom stellar population synthesis technique whose 
description is deferred to Section [3j Longslit spectro- 
scop ic d ata from numerous sources are described in Sec- 
tion [2?3 along with the derivation of gas-phase metallic- 



ities and Ha-based star formation rates. 

2.1. Supernova Sample 

The SNe la whose hosts are analyzed here were ob- 
served as part of the ongoing science operations f or the 
Nearby Supernova Factory (SNfactory - Aldering et al. 
2002[ ). The SNfactory was designed to observe sev- 
eral nundred SNe la in the nearby smooth Hubble flow 
(0.03 < z < 0.08) with the goals of achieving a deeper 
physical understanding of SNe la, building better SN la 
templates for cosmological applications, and anchoring 
the low-redshift Hubble Diagram. We briefly summarize 
the details of these operations. 

From 2005 to 2008, the SNfactory conducted a wide- 
field search of the northern and equatorial sk y using 
the QUEST-11 CCD camera ( [Baltay et aL]|2007[ ) on the 
Samuel Oschin 1.2m Schmidt telescope on Mount Palo- 
mar, California, partly in collaboration with the JPL 



Near-Earth Asteroid Tracking (NEAT) component of 
the Palomar-QUEST consortium. Typical search images 
consisted of 60 s exposure with an RG610 filter, with 
each field revisited multiple times in order to detect as- 
teroids (which we reject). The SNfactory search covered 
an average unique area of 600 deg^ per night and cov- 
ered nearly half the sky (w20,000 deg^) each year. In 28 
months of searching, the SNfactory discovered over 1000 
supernovae of all types, and spectroscopically confirmed 
over 600 of those. 

Spectroscopic typing of search candidates and followup 
observations of SNe la were obtained wit h the SuperNova 



Integral Field Spectrogr aph (SNIPS - Aldering et al. 



2002[ Lantz et al. 20041, mounted continuously on the 
University ot Hawan 2.2-m telescope on Mauna Kea. In 
addition to its spectroscopic capabilities, SNIFS has an 
imaging channel which was used to obtain broadband 
photometry of som e SN la host galaxies in our sample 
(see Section 2.2.1 1. 

A total of 400 supernovae discovered with the SNfac- 
tory search were spectroscopically confirmed to be 
SNe la, and those SNe la discovered before B-band max- 
imum light (as estimated by spectroscopic typing) were 
followed up extensively with SNIFS. In addition to those 
SNe la discovered by SNfactory, some SNe la discovered 
by other searches were followed with SNIFS, especially 
during times when search operations were impacted by 
weather or fires on Mount Palomar. This work analyzes 
all SNe la discovered or followed by SNfactory, a total of 
469 SNe la observed from 2004-2010. Of these SNe la, a 
total of 216 were extensively followed with SNIFS, with 
185 of those being discovered by the SNfactory search. 

In Table [T] we present the location of the SN la host 
galaxy as det ermined from the g-band host image (see 
Section [2.2.2[ ), the redshi ft measured from the host spec- 
trum (see Section 2.3.2) for all SNe la in our sample, 
and the discovery source for SNe la not discovered by 
SNfactory. Most of the SNe in our sample which were 
not discovered by SNfactory were discovered by ama- 
teur ast ronomers, the Palomar Transient Factory (PTF 
Rau et al. 2009 | , the Lick Obs ervatory Supernova Search 
ITUSS Fifippenko et al.|2001|), or the ROTSE superno va 
search ( | Akerldlt^r|pomnQuh^|2006l [Yuan|[20Tol ) . 

A total ot 40 SJNe m our sample have been presented 
in other SN la samples. 2 7 of our sample were also part 



of the CfA SN la sample (jHicken et al. 2009 2012), 24 



were observed by the C arnegie Supernova Proj ect (C S P 



Contreras et al. 2010 Stritzin ger et al.||2011 ), 21 were 
observed by LOSS ( Ganeshal ingam et al.t 2 010 |), and 1 



(SN 2005hc) was part ot the [Gupta et al.| ( |2011 ) sam- 
ple. Host galaxy masses were pre viously aeri ved for 1 
(SN 2006X) SN la host by JNeill eFall ([20091, but we 



did not derive a host mass for this galaxy due to shred- 
ding in photometric source ext raction. 14 of our SN la 

and 



(2010) 



hosts have masses derived by Kelly et al. 

we found our values (see Section [3]) had good agreement 

with theirs, having a mean offset of 0.04 dex in stellar 

mass. 

2.2. Host Galaxy Photometric Data 

Photometric data for SNfactory SN la host galax- 
ies was gathered from public sources as well as our 
own targeted observations. Optical photometry wa s col- 
lected from the Sloan Digital Sky Survey (SDSS, York 



et aL||2000| Eighth Data Release (DR8, [Aihara et ah 
2011| ). NIR i mages from the Two Micron All Sky Sur- 
vey (2MASS,'Skrutskie et al.|20"06 | were obtained at the 
NASA/IPAC Infrared Science Archive (IRSi^, and for 
select hosts NIR data were available f rom tHe UKIRT 
Infrar ed Deep Sky Survey (UKIDSS, [Lawrence et al. 
20071). UV dat a taken with the Gala xy Evolution Ex- 



plorer (GALEX, Morrissey et al.|2007l satellite were ob 



tained from the GALEX online data archive at MAStF^ 

The public photometric coverage of our hosts was very 
good. 344 of the 469 hosts (73%) fell within the SDSS 
photometric footprint, all 469 hosts had 2MASS imaging 
data, 58 (12%) hosts had UKIDSS data, and 408 (87%) 
have GALEX AIS (All-Sky Imaging Survey - a shallow 
imaging survey) images. Additionally, 102 (22%) of our 
hosts have deeper GALEX imaging, mostly from the MIS 
(Medium Imaging Survey). 

The typical photometric depth for these surveys (for 
this work, this limit is effectively where the flux errors 
reach 5-10% for point sources) are 20th magnitude for 
SDSS, 17th magnitude for 2MASS, 20th magnitude for 
UKIDSS, 19th magnitude for GALEX AIS, and 21st 
magnitude for GALEX MIS. After these cuts, the to- 
tal number of hosts with data from each survey is 323 
(69%) from SDSS, 394 (84%) from GALEX, 272 (58%) 
from 2MASS, and 49 (10%) from UKIDSS. 

For those hosts without optical photometry from 
SDSS, we used SNIFS in imaging mode to obtain optical 
images of the host after the SN had faded (at least one 
year later). SNIFS was also used to obtain deeper pho- 
tometry for those faint hosts whose SDSS images were 
not deep enough (typically for nig > 19.0). A total of 
283 hosts (60% of the total sample of 469 hosts) have 
SNIFS photometry. For 10 hosts, g-band photometry 
was obtained with Kec k LR IS prior to spectroscopic ob- 
servations (see Section 2.3) of the hosts, and was later 
zero-pointed to either SDSS or SNIFS photometry. 

Below we describe our reduction of the SNIFS pho- 
tometry and our method of combining multi-band imag- 
ing data to obtain accurate common aperture photom- 
etry. Our means of deriving galaxy stellar masses and 
star formation rates from photometry is described in Sec- 
tion [3| A mosaic image of several SNfactory host galax- 
ies is shown in Figure[l] utilizing both SDSS and SNIFS 
optical data and spanning a large range (4 orders of mag- 
nitude) of stellar masses. 

2.2.1. SNIFS Photometry 

For those hosts without publicly available optical pho- 
tometry from SDSS, or those faint hosts for which the 
photometric depth of SDSS was insufficient, we obtained 
optical photometry using SNIFS in imaging mode. The 
photometric imager ( "P-channel" ) on SNIFS consists of 
two 2kx4k E2V CCDs, with one "guider" chip undergo- 
ing fast continuous windowed readout to perform guiding 
during observations, and the other "science" chip dedi- 
cated to photometry. In normal SNIFS SN la observa- 
tion mode, the P-channel uses our custom "multi-filter," 
which observes field stars on the science chip segmented 
into multiple filters spanning a broad wavelength range, 
in order to estimate atmospheric extinction. The SNIFS 

^® http://irsa.ipac.caltech.edu 
^"^ http://galex.stsci.edu 



Table 1 

SNfactory SN la Host Galaxy Properties 



SN Name 



Host Coords 



Rcdshift (Hclio.) R.cfcrcncc Discoverer 



SNF20050519-000 


21:35:52.58 


-26:27:06.1 


0.03040 ± 0.00500 


SNF 


SNF20050621-001 


20:45:17.18 


-03:49:37.5 


0.12072 ± 0.00011 


SNF 


SNF20050624-000 


17:34:46.48 


-02:21:26.2 


0.06719 ±0.00009 


SNF 


SNF20050704-008 


22:26:29.49 


-16:41:38.4 


0.11116 ±0.00005 


SNF 


SNF20050727-005 


21:09:08.78 


-15:38:45.9 


0.08610 ± 0.00007 


SNF 


SNF20050728-000 


00:57:29.59 


-22:35:15.3 


0.04293 ± 0.00009 


SNF 


SNF20050728-001 


00:11:17.78 


-28:54:45.3 


0.05651 ± 0.00010 1 


SNF 


SNF20050728-006 


21:58:04.87 


-14:00:01.8 


0.05921 ± 0.00007 


SNF 


SNF20050728-012 


23:56:29.32 


-22:49:32.1 


0.07477 ±0.00250 


SNF 


SNF20050729-008 


22:06:56.95 


-09:29:30.4 


0.08206 ±0.00005 


SNF 



References: (1) Smith et al. ( 2004 1 
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Figure 1. A subsample of tlie SNfactory host galaxies, presented in gri color composites constructed using STIFF ( |Bertin|2012[ l. Galaxies 
are order by stellar mass from highest (upper left) to lowest (lower right). 



P-channel is also equipped with a variety of broadband 
filters covering the full science chip, including the stan- 
dard Gunn ugriz filters employed by SDSS. Because 
the SNIPS P-channel CCDs are different from those on 
the SDSS imager, the effective SNIPS filter throughputs 
vary slightly from those of SDSS. We show in Pigure [2] 
the fiducial SNIPS filter throughputs derived from the 
throughput of all the optical comp onents involve d and 
our Mauna Kea extinction curve ( Buton et al. 20136 com- 
pared to the SDSS filter throughputs. 

SNIPS images were processed in IRAlf^ using stan- 

^* IRAF is distributed by the National Optical Astronomy Ob- 



dard techniques. Overscan subtraction was performed 
for both amplifiers on the science chip, and data from 
each amp scaled by its gain. Images were trimmed to 
remove vignetting by the filter casing, then divided by 
normalized flatfield dome images to remove pixel vari- 
ations in detector efficiency. Por the reddest filters (i- 
and z-band), fringe patterns were removed by scaling a 
master fringe frame to the fringing measured in sky pix- 
els for each science image. Master fringe frames were 



servatory which is operated by the Association of Universities for 
Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 
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Figure 2. SNIPS filter throughputs, compared to those of SDSS. 

constructed from numerous long exposures, and identi- 
fication of sky pixels and fringe scaling were performed 
using custom soft ware. Cosmic rays were then removed 
using LA Cosmic ( van Dokkum|200T| . Astrometric solu- 
tions for all images were derived u sing WCSTo ols (MinP 
2006| , then refi ned using SCAMP (feertin| [2006J niatch- 
mg to 2MASS ( [Skrutskie et al.|[200S| ~astrometry. Im- 
ages from fie lds with multiple e xposures were combined 
with SWARP (Bertin et al. 2002) using median addition. 
For the purpose of deriving photom etric z eropoints and 
host galaxy photometry (see Section 2.2.2), fluxes for all 
objects in each image w ere measured with SExtractor 
( [Bertin fc Arnouts|1996l) using the FLUX_AUTD parameter, 
'i'he observing priorities for the SNIFS host photome- 
try program were to obtain g-band and z-band photome- 
try of all our hosts. The optical g — i color is a very good 
color for determin ing galaxy mass-to-light ratios (|Gal- 



lazzi fc Bell 20091 and thus serves as a minimal filter 
set tor obtainmg accurate galaxy masses. Many obser- 
vations were taken between the Seventh Data Release 
(DR7, Abazajian et al. 2009) and Eighth Data Release 
(DR8, jAihara et af:pnTp~ bf SDSS, which added a sig- 
nificant area to the bUbS imaging footprint. Thus we 
have a large number of fields in the SDSS footprint ob- 
served by SNIFS, especially in g and i, and with many of 
those observed on photometric nights when photometric 
calibration solutions were derived. This enables both the 
study of SNIFS-SDSS color terms as well as an indepen- 
dent measurement of the accuracy of our photometric 
calibrations, and we describe these two studies below. 

Photometric zeropoints for SNIFS imaging in the SDSS 
footprint were obtained by matching photometric mea- 
surements of field stars from each SNIFS science im- 
age to their values in the SDSS DR8 catalog. Formal 
zeropoints and their uncertainties were derived as the 
weighted mean (weighted by photometric error) of the 
zeropoints for individual field stars after the exclusion 
of severe outliers. In Table [2] we summarize the total 
number of SDSS fields visited with SNIFS in each band 
N fields, the mean of the zeropoint uncertainties {(Tzp) for 
all SNIFS fields zeropointed with SDSS in that filter, and 
the total number of stars matched over all fields Nstars ■ 
Most of our zeropoints are derived from w 30 field stars 
per image and have a precision of 0.01-0.02 mag. 

As stated above, the filter throughputs from SNIFS dif- 



fer slightly from those of SDSS, so we might expect small 
but nonzero color terms between the two filter sets. We 
can measure these from the same field stars used for ze- 
ropointing SNIFS photometry in the SDSS footprint. To 
do so, we compare the residual magnitude offsets (after 
application of the fitted zeropoint) of these field stars as 
a function of their color as measured by SDSS. We de- 
rive the weighted mean offsets in bins of color (typically 
0.2 mag wide) and perform a minimization to derive the 
optimal color term and its uncertainty for each filter. 
These are summarized in Table [2] As can be seen, the 
color terms are consistent with zero for all of the filters 
except z-band, which has a small but significant detection 
of a color term. This may be due to different strength 
of telluric absorption features at the SNIFS site (Mauna 
Kea) compared to the SDSS site (Apache Peak), or may 
be due to the slightly different red wavelength roll-off of 
the filter throughputs, as both of these effects may alter 
the shape of the filter throughput and drive a color term 
between SNIFS and SDSS. 









Table 2 










SNIFS color terms 




Filter 


N fields 


Nstara 


(o-zP> 


Color Term 


Color 


u 


9 


192 


0.0185 


-0.0009 ±0.0269 


u-g 


9 


160 


4914 


0.0094 


0.0004 ± 0.0087 


g-r 


r 


12 


790 


0.0109 


0.0014 ±0.0104 


g-r 


i 


157 


12452 


0.0143 


-0.0222 ±0.0115 


r — i 


z 


12 


1068 


0.0294 


0.0081 ± 0.0561 


i — z 


9 


- 


- 


- 


-0.0010 ±0.0044 


9-i 


i 


- 


- 


- 


0.0099 ± 0.0056 


9-i 



Photometric zeropoints for SNIFS fields outside the 
SDSS footprint were derived for each observing night in 
each filter using observations of standard stars spanning 
a large r ange of airm asses. Our standards were selected 
from the Smith et al.| ( |2002[ ) sample, placing our measure- 
ment on the standard ugriz system employed by SDSS. 
For each night (in each filter) we fit for a global zeropoint 
and an atmospheric extinction term, and our extinction 
terms were consistent with tho se predicted by the fiducial 



Mauna Kea extinction curve (Buton et al. 2013). Typi 



cal dispersion of standard star inagnitudes about the best 
fit calibration solution were about 0.02 mag in gri and 
0.03 mag in u and z. New science images were assigned 
a zeropoint based on their airmass and exposure time as 
calculated with the fitted extinction solutions. 

As stated above, a number of the fields for which we 
obtained new zeropoints were included in the subsequent 
SDSS data release, enabling us to derive external zero- 
points to cross-check our calibration solutions. We com- 
pared the SNIFS-based zeropoints to those derived by 
matching to SDSS and found good agreement (mean ze- 
ropoint offsets less than about 0.005 mag) with a disper- 
sion consistent with the dispersion seen in our calibra- 
tion solutions (about 0.02-0.03 mag). Since the SDSS 
zeropoints are derived from a large number of field stars, 
they are ultimately more precise than those derived from 
SNIFS observations of a few standard stars. We thus 
use SDSS zeropoints in favor of SNIFS zeropoints where 
available. 

2.2.2. Final Host Galaxy Photometry 
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Figure 3. Example of common aperture photometry for the host 
of SNF20060609-002, showing the resampled images in each band 
and the photometric aperture. 



With the final processed SNIPS imaging and pubUc 
data from SDSS, 2MASS, UKIDSS, and GALEX, we ob- 
tain magnitudes for our hosts in each band by performing 
common aperture photometry. We use the f;-band im- 
age to select the host galaxy and define the galaxy aper- 
ture. For each SN, the location and elliptical aper tures of 
nearby galaxies were measured with SExtractor (Bertin 
|& Arnouts 1996) and the separation of each galaxy from 
the SN location was calculated and then scaled to the ef- 
fective galaxy radius projected along the SN-galaxy vec- 
tor. The SN host was then nominally identified as the 
galaxy the least number of effective radii from the SN 
location. A spectroscopic redshift of the prime galaxy 
candidate was required to match the SN rcidshift as de- 
termined using SNID (Blondin fc Tonry||2007j ) to within 
1500 km s~^. In nearly all instances the nearest host 
candidate was confirmed to be the host, though there 
were a few cases in which the nearest host candidate was 
a background high-redshift galaxy. 

For some of the PTF targets from 2010, wc will re- 
port preliminary host properties in Section [3] based on 
the SN redshift and the nominal host galaxy, some SNe 
appear to be hostless, and these interesting cases will be 
a subject of future study. 

After positive host identification, we measure the host 
flux in each band by first resampling the image from each 
filter to the resolution of the ap erture image us ing SWARP 
and then running SExtractor (Bertin & Arnouts 1996| 
in dual image mode. We use the SExtractor FLUX_AUTO 
output parameter, which measures the flux inside an el- 
liptical Kron-like aperture. Sky background subtraction 
was performed using a large background mesh (1024 pix- 
els, equivalent to about 4 arcmin). Star contamination of 
galaxy photometry was handled with SExtractor object 
deblending, which assigns pixels to overlapping objects 
based on where the flux of the star drops below a cer- 
tain contrast threshold, which we chose to be the default 
fraction of 0.005. An example of our common aperture 
photometry method is shown in Figure [3] 

Galaxy magnitudes and their errors in each band were 
calculated using the zeropoints and noise characteris- 
tics for each image. We convert all magnitudes to the 
AB systems by applying Veg a-AB offsets fo r 2MASS 



and UKIDSS magnitudes (Hewett et al 



2006), and the 



ugriz offsets for SDSS as derived by Kessier et al.| ( |2009[ ) 
GALEX zeropoints are already on the AB system. Ub- 
served magnitudes were then corrected for fo reground 
Milky Way reddening using the dust m aps of Schlegel 
et al. (1998) and the reddening law of Cardelli et al. 



( [19891 hereafter CCM). 

Fot some images, such as GALEX UV images of ellip- 
tical hosts with no star-formation, the galaxy flux was 
below the image noise threshold. For these images we 
calculate an effective 3cr magnitude upper limit deter- 
mined by the flux uncertainty measured in the galaxy 
aperture after accounting for Milky Way dust reddening. 
These upper limits, especially those from GALEX UV 
data, proved valuable in constraining the stell ar p opula- 
tion flts to photometry, as outlined in Section [3. 4| 

Finally, optical magnitudes were corrected for emission 
line fluxes using the fitted line profiles and stellar con- 
tinuum mea surements obtained from optical spectra (see 
Section 2.3) where available. These emission corrections 



were typically less than 0.01 mag in all bands except 
for r or i (depending on redshift) where the strongest 
emission line (Ha) was typically present. For the most 
strongly star-forming hosts in our sample, the emission 
corrections still were less than 0.05 mag. 

To convert these galaxy photometric data into physical 
properties of the galaxies, one must invoke the art of 
stellar population synthesis (SPS). Because this is a rich 
and constantly evolving field, we devote the entirety of 
Section [3] to describing our chosen method of deriving 
galaxy stellar masses and SPRs. 

2.3. Host Galaxy Spectroscopic Data 

Galaxy spectroscopy is useful for gaining finer insight 
into the galaxy SED than can be gleaned from broad- 
band photometry. In particular, absorption features in 
the stellar continuum of the galaxy SED can be com- 
pared to stellar evolution models to estimate stellar age 
and metallicity, while narrow emission lines from ionized 
HII regions surrounding young stars can yield both gas- 
phase metallicity of the galaxy interstellar medium (ISM) 
as well as the current rate of star formation. Addition- 
ally, reddening in the galaxy ISM and near the ionized 
HII regions can be estimated from galaxy spectra. In 
this Section we describe both the SNfactory host galaxy 
spectroscopic data set as well as the extraction of galaxy 
physical parameters from these data. 

2.3.1. SNfactory Host Spectroscopy Observations 

Longslit spectra for our SN la host galaxies were ob- 
tained during numerous observing runs at multiple tele- 
scopes from 2007-2011. The instruments used were the 



Kast Double Spectrograph (Miller fc Stone|1993) on the 



Shane 3-m telescope at Lick Ubservato ry, the Low Res- 



Oke et al. 



1995 1 
the R- 



olution Imaging Spectrometer (LRIS 
on the Keck I 10-m telescope on Mauna Kea, 
C Spectrograph on the Blanco 4-m telescope at Cerro 
Tololo Inter-American Obs ervatory, the G oodman High 
Throughput Spectrograph ( jClemens et al.„2004J on the 
Southern Astrophysical Research ( SOAR) 4-m telesc ope 
on Cerro Pachon, and GMOS-S (jDavies et al.||1997[ ) on 
the Gemini-S 8-m telescope on Oerro Pachon. The in- 
strument configurations, including wavelength coverage 
and effective resolution (FWHM), are presented in Ta- 
ble [3j Objects were generally aligned onto the slit using 
direct imaging techniques where available, or offset star 
alignment when direct imaging was not available. The 
slit position angle was generally chosen to correspond to 
the parallactic angle at the time of observation, though 



occasionally when the airniass was favorable (i.e. less 
than about 1.05) the slit was aligned along the galaxy-SN 
direction to enable future abundance gradient or rotation 
studies. 

Longslit spectra were reduced in IRAF using standard 
techniques. After overscan subtraction, we subtracted 
bias frames from two-dimensio nal longslit data, rem oved 
cosmic rays using LA Cosmic (van Dokkum 2001), and 
flatfielded to remove pixel variations in detector effi- 
ciency. Two-dimensional wavelength solutions were de- 
rived from arc lamp exposures taken either at the same 
pointing as the object spectrum (for Shane, Blanco, and 
SOAR data) or using nightly arc lamp exposures (for 
Keck and Gemini-S data), with a one-dimensional shift 
applied by measuring atomic (01) night sky lines in ob- 
ject spectra. Object spectra were reduced to one dimen- 
sion using the IRAF function apall, with the galaxy 
aperture chosen by visual inspection to encompass as 
much of the galaxy core as was available above the sky 
noise. Nightly flux calibrations were derived from stan- 
dard stars observed at appropriate ranges of airmass, 
with an atmospheric extinction solution derived for each 
observing night. Telluric absorption features were then 
removed using the nightly standard star spectra, ac- 
counting for appropriate dependence on airmass. Ob- 
server motion with respect to the heliocentric frame was 
then corrected, and finally spectra were dereddened to 
correct for Milky Wa y extinction using the dust maps of 



Schlcgcl e t_al.| ( 1998 ) and the reddening law of Cardelli 
et ai.,(.198 



Some hosts had spectra available from SDSS DR8 ( Ai 



hara et al. 20111. These spectra were downloaded and 



then converted to air wavelengths for consistency with 
reduction of our own observations. The total number of 
spectra from each source was 225 host spectra from Kast, 
82 from SDSS, 29 from the R-C Spectrograph, 18 from 
LRIS, 13 from the Goodman HTS, and 7 from GMOS- 
S, for a total of 374 host spectra (81% of the sample 
with host photometry) . In Table [4] we list the instru- 
ment, UTC date of observation, exposure time, slit PA, 
and mean airmass for each SN la host spectroscopic ob- 
servation. For SDSS spectra, we list the plate-mjd-fiber 
combination of the spectroscopic observation in the ob- 
servation date column. 

We note here that our visual selection of spectrum ex- 
traction apertures may result in a different distribution 
of spectroscopic covering fraction as compared to other 
spectroscopic surveys such as SDSS. For example, the 
typical extraction aperture used for our Kast observa- 
tions was 10" with a 2"-wide slit, whereas the SDSS spec- 
troscopic fibers covered a 3" diameter circular aperture. 
For the metallicity and SFR comparisons undertaken in 
Sections |4] and [5] we are interested in how our aperture 
selections might affect the metallicity and SFR values 
measured from emission line fluxes in our host galaxy 
spectra. 

Differences between galaxy metallicity or line emission 
measured with different ape rtures have been stu di ed by 
a number of authors, e.g., Gom ez et aT] (|2003|); Kew- 

(|2007); Moustakas et al. 
nese differences depend on 
covermg traction, but m different ways for different mea- 
surements. In particular, Ha flux can vary strongly with 
position due to HII regions, thus affecting Ha-based SFR 



ley et al.| (2005); Salim^et aT 
JU10| ; |Zaliidet ai.| ( ,20f^ rT 
lovermg fraction, but in differ 



and equivalent width (EW) measurements, while metal- 
licities will vary less strongly with positio n, e.g. primar- 
ily du e to radial gradients. For instance Gerssen et al. 
(20121 find that SDSS near-core Ha emission has a poor 
correlation with global Ha for captured-light fractions 
below 40%. [Kewley et al.| p005) find that for slit spec- 
troscopy a captured-light fraction greater than ^20% is 
adequate to achieve agreement between global and slit 
metallicities. 

For a more quantitative estimate of the effect that 
galaxy metallicity gradients and aperture choice may 
have on our re s ults, we turn to the recent study by 



Sanchez et al. (2012|. They found that metallicity 
gradients in spiral galaxies are very consistent across 
all spiral galaxy sub-types when scaled by the effec- 
tive radius (rg) of the galaxy, with mean and scatter of 
—0.12 ± 0.11 dcx/rg. Our apertures cover typically one 
scale radius, so the light-weighted metallicity for an ex- 
ponential disk would be only 0.03 ± 0.03 dex lower than 
the at-core metallicity. Since our metallicities are mea- 
sured from emission lines rather than starlight, the light 
profile of the e mission provides a more accurate result. 
Koopmann et al. ( 2006 ) found the Ha scale radius to be 
1.14X the stellar scale radius on average. This implies 
that our emission-weighted metallicities would be only 
0.04 ± 0.04 dex lower than the core metallicities. Even 
if integrated to infinite radius, our metallicities would 
be only 0.08 ± 0.08 dex lower than core for a emission- 
weighted extraction of an exponential disk. These poten- 
tial offsets are well below the systematic uncertainties on 
the metallicity methods themselves and thus are unlikely 
to influence our results. This notion is supported by the 
fact that our hosts show such excellent agreement with 
the galaxy mass-metallicity relation (Section 111). 

2.3.2. Redshifts and Emission Line Fluxes 

SN la host galaxy redshifts, metallicities. Ha star- 
formation rates, and internal reddening were calculated 
using emission line fluxes from the host galaxy spectra. 
Accurate measurement of emission line fluxes in star- 
forming galaxies requires proper accounting for stellar 
absorption. To this end we fit the emission line fluxes and 
stellar background in each host spectrum simultaneously 
using a modified version of the IDL routine linebackf it 
from the idlspec2q^ package developed by the SDSS 
team at Princeton. This routine allows the user to pro- 
vide a list of template spectra which are then fitted to 
the data in linear combination with Gaussian emission 
line profiles. We have modified this code to force the 
coefficients multiplying the stellar continuum templates 
to be non-negative and to fit for reddening of the stellar 
continuum (using a CCM law with Ry = 3.1 fixed and 
E{B — y) as a fit parameter). Additionally, we have in- 
corporated the ability to fit for a scaling factor between 
the blue and red channels of two-arm spectrograph data. 
For background templates we chose a set of simple stellar 
populations (SSP s) from the stellar popula tion synthesis 
code GAL AXEV (B ruzual fc Charlot|2 003', BC03) with a 
Chabrier ( 2003) IMI*' and the same age sampl ing used for 
background htting by Tremonti et al. (2004), which ul- 
timately consists of ten bbPs tor each metallicity track. 
These templates are convolved to the resolution of the 

^^ http://spectro.princeton.edu/idlspec2d_install.html 



Table 3 

Instrument Configurations 



Instrument 


Dichroic/ 


Disperser 


Slit 


Wavelength 


Effective 




Filter 




(arcsec) 


Coverage (A) 


Resolution (A) 


Kast blue 


d55 


600/4310 


2.0 


3900-5550 


3.1 


Kast red 


d55 


300/7500 


2.0 


5450-10500 


9.1 


LRIS blue 


D560 


600/4000 


1.0 


3500-5600 


3.9 


LRIS red 


D560 


900/5500 


1.0 


5500-7650 


4.2 


Goodman HTS 


GG385 


300 1/mm 


1.0 


3850-7700 


13.7 


R-C Spec 


GG385 


300/7500 


1.0 


3850-7700 


9.1 


GMOS-S 


GG455 


B600 


1.5 


5040-7920 


6.8 



Table 4 

SNfactory SN la Host Spectroscopic Observations 



SN Name 



Instrument 



UTC Date 



Exp. Time Slit PA Airmass 



SNF20050621-001 


Kast 


2007-Oct-15.2 


1800 


111 


1.499 


SNF20050624-000 


Kast 


2007-Jun-15.4 


1800 


82 


1.326 


SNF20050704-008 


R-C Spec. 


2009-Aug-28.2 


900 


100 


1.029 


SNF20050727-005 


R-C Spec. 


2009-Aug-30.1 


1200 


48 


1.091 


SNF20050728-000 


Kast 


2007-Oct-14.3 


1800 


178 


1.992 


SNF20050728-006 


Kast 


2010-Sep-11.3 


1200 


188 


1.614 


SNF20050728-012 


R-C Spec. 


2009-Aug-29.2 


900 


189 


1.010 


SNF20050729-008 


SDSS 


0718-52206-0264 


9211 




1.443 


SNF20050730-003 


Kast 


2009-Oct-24.2 


600 


182 


1.309 


SNF20050731-005 


Kast 


2009-Aug-21.3 


3600 


107 


1.314 



particular spectrogr aph w hose data we fit. We note that 
the use of Salpcter, ( 1955 ) IMF templates results in neg- 
ligible differences to the fitted emission line fluxes, and 
metallicity difference smaller than our typically quoted 
precision of 0.01 dex. An example fit to spectroscopic 
data is shown in Figure |4J 

Use of the continuum model naturally accounts for 
the effect of Balmer absorption. Given our typically 
low spectral resolution (e.g. 400 km s^^ at Ha and 
200 km s^^ at H/3 for Kast), we did not fit for a veloc- 
ity dispersion for the stellar continuum in each SN la 
host galaxy spectrum since our templates were already 
broadened by the spectral resolution. This is adequate 
for nearly all host galaxy and instrument combinations, 
where the expected galaxy velocity dispersion (based on 
Gy oc M^/'^) is either below the spectrograph resolution 
or Ha is strong enough that stellar absorption correc- 
tions are very small relative to the emission strength. 
For example, we re-fit the Kast spectrum of the strongly 
star-forming {EW{Ha) = 75A) host of SNF20050826- 
004 with stellar templates broadened by 400 km s~^, 
and found the Ha flux changed by less than 2%. There 
may be some massive ellipticals with weak emission (e.g. 
log(M*/M0) > 11 and EW{Ha) < 5A) for which veloc- 
ity dispersion could be important, but these cases typi- 
cally have emission dominated by AGN activity and ul- 
timately do not produce a gas-phase metallicity. 

To ensure that our modifications of the code are not 
producing anomalous results, we compared our fitted 
emission line fluxes to those derived by the MPA-JHU 
team for those hosts whose spectra were obtained from 
SDSS. In Figure[5] we plot these values and show that our 



results are very consistent with those derived by these au- 
thors. Thus we believe our emission line flux estimates, 
including the resultant corrections for continuum absorp- 
tion, are accurate. 

The only minor discrepancy may be in the flux of the 
H/3 line. Our emission line fits use the BC03 models, 
whereas the MPA-JHU values for DR7 use the updated 
(unpublished) Chariot & Bruzual 2008 models. The 
BC03 models have been previously shown to underesti- 
mate the strength of the H/3 emission line ( Groves et al. 



2012 1 in a comparison between the DR7 and UK4 fits 
(which used the BC03 models) due to inaccurate fitting 
of the H/3 absorption profile. This issue could affect esti- 
mates of the galaxy reddening and thus bias metallicity 
measurements. For example, an underestimate of the H/3 
line would over-estimate the reddening, which would re- 
sult in an over-estimate of the OH A3727 flux, which 
in the R23 metric would produce an over-estimate of 
the metallicity. However, we use a reddenin g-insen sitive 
metallicity indicator (PP04 A^2 - see Section 2.3.3) in all 
but the lowest metallicity SN la hosts. The H/3 bias (of 
order a few A in equivalent width) is strongest in galax- 
ies with weak emission, and our low-metallicity hosts are 
all very strongly star-forming (EW(H/3)~50A) and thus 
likely to be only marginally affected by this problem. 

Redshifts for SNfactory host galaxies with strong emis- 
sion lines were derived as the weighted (by measurement 
uncertainties) mean of individual emission line redshifts 
fitted from host spectra. Redshift errors were similarly 
calculated from the measurement uncertainties on the in- 
dividual line redshifts. This method is the same as that 
used by SDSS for final redshifts reported in their online 
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Figure 4. An example of our fits to galaxy spectroscopy, here for ttie liost of SNF20070331-013 which was observed with LRIS on Keclc. 
The blue curve is the data, the green is the stellar continuum fit, and the red is the fitted emission line profiles. In the lower panel we show 
a zoomed-in view of each of the fitted lines, along with the residual difference (magenta curve) between the data and fitted model of the 
stellar continuum plus emission lines. 

der the assumption that nebular emission is excited by 
ionizing flux from young stars. To do so we must first ex- 
cise those galaxies whose emission line fluxes are contam- 
inated by AGN activity us i ng the emission line diagnostic 
diagram of Baldwin et al. ( 1981 hereafter BPT). In Fig- 
ure^ we plot the distribiition of SNfactory emission line 
host galaxies on the BPT diagram as compared to the 
distributio n of galaxies from SPS S, with the boundaries 
defined by Kewley et al. (2006) to distinguish normal 
star-forming galaxies frorn AGJNs and composite galax- 
ies. Of the 374 host galaxies with spectra, 315 had good 
detections of all four lines used in this BPT diagram. Of 
those, 215 had emission line fluxes consistent with ion- 
ization by young stars, 50 had line fluxes indicating AGN 
ionization, and 50 had line flux ratios in the composite 
region of the diagram. This high AGN fraction is likely 
a product of our spectra being dominated by light from 
the galaxy core where AGN activity is strongest. 

Finally, for those hosts whose emission was consistent 
with star- format ion, emission line fluxes were corrected 
for internal reddening within the host galaxy by em- 
ploying the Balmer decrement method. In an HII re- 
gion ionized by young stars, the ratio of emission line 
flux in the Ha line to that in the II/3 line is fixed by 
atomic physics. For typical HII regions whose tempera- 
tures are around 10,000 K, this ratio is calculated t o be 
F(Ha)/ F(Hl3) = 2.87 under Case B recombination ([Os" 




Figure 5. A comparison of the fitted emission line fiuxes from our 
modified emission line fitting code (y-axis) vs. the values derived 
by the IMPA-JHU team (x-axis). Points are color coded by emission 
line, and the line representing unity is the solid black line. 

database. For hosts with very weak or no emission lines, 
redshifts were calculated with a cross- correlation tech- 
nique using the methods presented by [Tonry fc Davis] 
( 1979 1 . We correlated the best fit stellar continuum spec- 



trum against the observed host spectrum after subtrac- 
tion of the fitted emission line fluxes. Typical redshift er- 
rors for these two methods are of the order Cz ^ 0.0001, 
and tend to be limited by the precision of the spectro- 
scopic wavelength solution. 

In order to derive gas-phase metallicities (Sec- 
tion 2.3.3 ) or Ha-based star- format ion rates (Sec- 
tion 2.3.4), we will invoke empirical formulae derived un- 



terbrock fc Ferland||2006[ ), and is known as the Balmer 
decrement. Reddening by dust causes the observed value 
of this flux ratio to exceed its canonical value, and one 
can calculate the amount of reddening by assuming a 
reddening law such as that of Cardelli et al. (1989|. For 
our SN la host galaxies, final emission line fluxes used for 
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Figure 6. SN la hosts from SNfactory in the BPT diagram. The 
diffuse grey background represents SDSS galaxies whose emission 
line fluxes were measured by the MPA-JHU team. Galaxies above 
the solid red line are classified as "AGN" galaxies, those below the 
dashed blue line are "star-forming" galaxies free of AGN contam- 
ination, and those galaxies between the two lines arc classified as 
"composite" galaxies. 



calculations of gas-phase metallicity and star-formation 
rate (from Ha) have been corrected for the internal red- 
dening calculated using this method. 

2.3.3. Host Gas Phase Metallicities 

Translating emission line fluxes into a gas-phase metal- 
licity depends on the ch oice of metallicity ca,librat ion, as 
thoroughly detailed by [Kewley fc Ellison| ( |2008[ ). Dif- 
ferent calibrations are known to disagree by as much as 
~ 0.5 dex, which makes it difficult to place metallicity 
measurements on a common absolute scale. Addition- 
ally, there is no single metallicity metric that is ideal 
across the entire range of metallicity probed by our sam- 
ple. For exa mple, metrics that r ely on the Nil A6584 
line, such a s Kewley & Dopita| (|2002|) and |Pettini &| 
|Pagel (20041) methods, have high signal-to-noisc at high 
metallicity and are monotonic, but at low metallicities 
this line becomes very weak and produces large errors 
in metallicity measurements. The well-known R23 met- 



( 2004 herea fter T04) scale using the c onversion formulae 
presented in Kewley & Ellison ( 2008 1 . 

For galaxies with log(JMl/Haj > —1.3 (i.e. " high" 
metallicity ga l axies) , we use the "N2" method of Pet- 
tini & Pagel (2004), as the NII/Ha ratio is a sensi 



tive metaUicity indicator in this range and has rela- 
tively low sensitivity to reddening due to its short wave- 
length baseline. For very low metallicity galaxi es with 
log(NII/Ha ) < -1.3, we use the "R23" method of|Kobul 



nicky fc Kewley (2004) as updated by (Kewley & Ellison 
20081), as this rnethod depends on the relatively strong 



oxygen lines and also fits iteratively for the ionization 
parameter. Although the R23 metric is doubly valued 
with metallicity, the choice of our log(NII/Ha) cut places 
these galaxies firmly on the low-metallicity "branch." 
Metallicities calculated from these original methods are 
finally converted to the T04 scale u sing the conversion 
formulae of Kewley fc Ellison (2008). It is worth noting 
that the dispersion in these conversion formulae is com- 
parable to the systematic uncertainty in the metallicity 
calibrations themselves. Thus unless otherwise noted, all 
metallicities reported here are on the T04 scale after ap- 
plication of the above described conversion. We report 
the metallicity method and final T04 metallicity for all 
our SN la hosts in Table [5] along with the BPT clas- 
sification of our host galaxies, and the reddening value 
calculated from the Balmer decrement method. 

2.3.4. Host Star- Formation Rates from Ha Flux 

The flux in the Ha line is often used to calculate the 
current star- format ion rate (SFR) in galaxies. We calcu- 
late Ha-based SFRs for our star-forming (i.e. non-AGN) 
hosts using the following method. We first measure the 
Ha fiux from the spectrum, then correct it for redden- 
ing using the Balmer decrement as described above. To 
determine the total amount of Ha flux corrected for spec- 
troscopic slit-loss, we calculate a synthetic i-band mag- 
nitude from our galaxy spectrum and compare it to the 
observed i-band magnitude to calculate a scaling factor. 
From this total observed galaxy Ha flux, we calculate the 
total rest-frame Ha luminosity using the distance niod- 
ulus at t he ga laxy's redshift calculated with the code of 



ric (e.g. McGaugh 1991 Zaritsky et al. 1994; Kobulnicky multiply by 0.7 to convert t he Kennicutt for mula 



k, Kewley|2004) is double- valued with metallicity, and is 



Wright] ([2006 ). Finally we con vert this Ha lumi nosity 
to a SFR using the formula of Kennicutt (1998) then 



rom a 



shallow-sloped (in terms of R23 versus metallicity) at low 
metallicity (i.e. flux errors propagate into larger metallic- 
ity errors). At low metallicities, the preferred me tallicity 



1984), which 



calibration is the re(OIII) method (Aller 
relies on the auroral A4363 oxygen line. 'Ihis method is 
considered the most reliable, but relies on a very weak 
emission line and does not consistently agree with the 
empirical strong-line methods. 

Thus it is challenging to flnd a consistent metallicity 
calibration that has high sensitivity over the full observed 
range of galaxy gas-phase metallicities. Additionally, 
the lack of consistency of absolute metallicity scales be- 
tween various calibrations makes it difficult to compare 
reported metallicity values from numerous authors. In 
order to utilize the strongest available lines and place 
our measurements on a well-known common scale, we 
employ different calibrations at different scales a nd then 
place all our metallicities on the common |Tremonti et al 



Salpeter] ([T955| IMF to ou r jChabrierl ( J2003l ) IMF. 

We can compare these Ha-based SFRs to those cal- 
culated from photometry using our SPS methods (see 
Section p|, and plot this comparison in Figure [7] In the 
same figure, we plot the SFRs for SDSS galaxies cal- 
culated from Ha fiux measured by the MPA-JHU team 
coupled to the same |Kennicu "tt] (jl998| formula against 
SFRs calcul ated by that gro up from broadband galaxy 
photometry ( Brinchmann et al.,2004) . As can readily be 
seen, the relation between spectroscopic Ha SFRs and 
photometric SFRs is very tight for SDSS, with a dis- 
persion about the unity line of 0.22 dex. For SNfactory 
hosts the agreement is also good but with an increased 
scatter of 0.66 dex which is larger than the typical mea- 
surement error of 0.23 dex for the Ha-based SFR values. 
The smaller scatter in the SDSS data is likely due to the 
fact that the MPA-JHU SFR values are the sum of the 
Ha-based SFR for galaxy fiux in the SDSS spectroscopic 
fiber and a color-based SFR estimate for galaxy fiux out- 
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Table 5 

SNfactory SN la Host Spectroscopic Results 



SN Name 



BPT Class Z Method 



Zt04 



+0.04 



Balmer E{B - V) 


EW(Hq) 


log{F(Ha)) 


0.058 ±0.145 


19.51 


41.34 




30.00 


41.65 


0.557 ±0.082 


22.47 


41.28 




10.53 


41.00 


0.000 ± 0.044 


47.27 


40.50 




10.91 


40.53 


0.933 ±0.151 


19.63 


41.07 




1.47 


40.32 




6.12 


41.02 




0.97 


40.41 



SNF20050621-001 


SF 


N2 {PP04) 


SNF20050624-000 


AGN 




SNF20050704-008 


SF 


N2 (PP04) 


SNF20050727-005 


AGN 




SNF20050728-000 


SF 


N2 {PP04) 


SNF20050728-006 


Comp. 




SNF20050729-008 


SF 


N2 {PP04) 


SNF20050730-003 


AGN 




SNF20050731-011 


Comp. 




SNF20050820-004 


Comp. 





9.11 



9.02 



-0.03 
-0.03 



8.36 



+0.08 



;.98 



+0.03 




-2-10 1 

SFR from Photometry 



Figure 7. Star-formation rates for SNfactory SN la hosts (blue 
circles) calculated from Ha flux in spectra versus those calculated 
from broadband photometry ( Section [Sjl. The same quantities cal- 
culated for normal star-forming galaxies from SDSS are shown as 
the grey density plot, while a line of unity is shown as the thick 
red line. 

side the fiber trained on fiber colors and Ha fluxes for the 
full spectroscopic sample. A more detailed discussion of 
these effects will be given in Section [5] We note that the 
SFRs of the SNfactory sample extend much lower than 
the nominal SDSS sample plotted here, due to the greater 
representation of low mass galaxies in the SNfactory host 
sample. 

3. HOST GALAXY MASSES AND 
STAR-FORMATION RATES FROM 
PHOTOMETRY 

Calculation of a galaxy's stellar mass and star forma- 
tion rate from its photometry requires the use of stellar 
population synthesis (SPS) techniques. The core princi- 
ple of this technique involves using model stellar popula- 
tion spectral energy distributions (SEDs) to predict the 
flux in various photometric filters, then comparing these 
model predictions to observations. SPS techniques typi- 
cally combine model SEDs for stars of a single age with 
masses distributed according to some initial mass func- 
tion (IMF), thereby deriving the SED for what is known 
as a simple stellar population (SSP) of uniform age and 
metallicity. Full galaxy SEDs are calculated by prepar- 
ing a model star- format ion history (SFH) and convolving 
the SSP SEDs with the relative weights prescribed by the 



galaxy SFH. 

The field of galaxy stellar population synthesis is a 
rich and constantly evolving field. The best SPS mod- 
els require stellar evolutionary tracks as well as observed 
(and/or modeled) stellar SEDs spanning the full param- 
eter space of stellar properties. While most SPS tech- 
niques give qualitatively similar results, it is important to 
understand and track the differences between SPS tech- 
niques employed by different authors. Perhaps the two 
most popular sets of models in the past decade have been 
GALAXEV dB ruzual fc C harlot ||2003l hereafter BC03) 
and PEGASE (I'ioc & Rocca-Volmerange||1997|). Most 



of the major analyses of the SDSS galaxy sample em- 
ployed the BC03 models (iKau ffmann ct al. 2003; Brincl>] 
' mann et al |2004||Tremontret al.|2004j | Sahm et al. | 2 007|7 



Many SIN la host galaxy studies from recent years ([SuiiT 



|2004j | 
from r 



van et al. 2006; HoweU ct al. 2009; Ncill et al. 2009; Kelly 
et al.|2010 ; Sullivan ct al..2010 , Lampeitl et al. 2010 
made use of PEGASE, but fortunately the diiteren 



have 
Breht SPS 

models gi ve consistent resu lts when scaled appropriately 
(see, e.g., Kelly et al.|2010 ). These particular SN la host 
studies employed the code ZPEG (Le Borgne & Rocca- 
|Volmerange ^2002) , which matches observed photometry 
to PEGASE models for a set of galaxy evolut i onary sce- 
narios input by the user [t hough Gupta et al. ( 2011 ) use 
the Flexible SPS code of [Conroy fc GunnJlMuf and 
a probabilistic galaxy parameter estimation approach in 
some ways similar to ours]. Though designed primarily 
as a tool for deriving photometric galaxy redshifts, ZPEG 
inherently derives galaxy masses and star- format ion rates 
by choosing a best model SED and scaling it to the ob- 
served galaxy photometry. 

For our purposes we desire not only an accurate esti- 
mate of our host galaxy properties (stellar mass, SFR, 
reddening) but also a proper estimation of how well the 
data constrain those properties (i.e. their uncertainties). 
This requires that we calculate not only how well a single 
model galaxy SED is fitted by our observed photometry, 
but also what range of models can fit our data and what 
this implies for the allowable range of galaxy properties. 
This inherently requires that we employ a Bayesian ap- 
proach, which we will describe in detail in the following 
section. 

3.1. Bayesian Stellar Population Synthesis of SN la 
Host Galaxy Photometry 

Deriving the physical properties of a galaxy from 
broadband photometry is a challenging task. Galaxies 
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are composed of billions of stars of various masses and 
ages, embedded within some nontrivial geometric distri- 
bution of gas and dust. A small discrete set of observ- 
able quantities, namely broadband photometry, give in- 
sufficient degrees of freedom to independently constrain 
the full stellar population and dust distribution within a 
galaxy. Instead we must establish some reasonable prior 
constraint on how the ages of stars within that galaxy are 
likely to be distributed by positing the functional form of 
the galaxy star-formation history (along with a similar 
prior on dust) , as well as the mass distribution produced 
by each generation of stars (the IMF). With this SFH 
prior and a chosen IMF we then use the observational 
data for a particular galaxy to infer the most likely star- 
formation history as well as its uncertainty. In practice, 
we implement this Bayesian inference technique by us- 
ing SFH and dust priors to produce a large sample of 
Monte-Carlo generated SFHs, which then define model 
SEDs when coupled to the SPS models of BC03. We then 
constrain three key properties of our observed galaxies - 
mass-to-light ratio, specific star- formation rate (sSFR), 
and dust reddening - by summing over the values of those 
properties for each model galaxy, weighting by how well 
the model photometry matches the observed photometry. 
Our method is inspired b y, and closely resembles , 
the methodolog y empl oy ed by Kauffmann et al. J 2003 (^ 



Gallazzi et al. 



([20051 ); ISalim et al.| ( |i>007| ; [Gallazzi 



^ Bel l (2d^:'^ Werandomiy generate LSTCoOO model 



galaxy Sl'Hs as follows. Model galaxy SFHs consist of 
an exponentially declining star-formation rate (SFR oc 
exp(— 7t)) with the decay constant (7) distributed uni- 
formly between and 1 Gyr"-'^ and star- formation begin- 
ning at a time iform distributed uniformly between 1.5 
and 13.5 Gyr in the past. In addition to the exponen- 
tial continuous star- format ion component, each SFH is 
augmented with random bursts of star formation whose 
durations are distributed uniformly in the range 30 Myr< 
iburst <300 Myr, with mass distributed logarithmically 
between 0.3 and 4.0 times the total mass formed in the 
continuous SF component. These bursts of star forma- 
tion occur with a random likelihood of occurrence such 
that 50% of model galaxies have experienced a burst of 
SF in the past 2 Gyr. Each model galaxy is assigned 
a single metallicity (i.e. no chemical evolution is mod- 
eled), with all model metallicities distributed uniformly 
between 0.2 < Z/Zq < 2.5 and as a smoothly decaying 
function (a log(Z)i/3) in the range 0.02 < Z/Zq < 0.2 
(in order to not over-represent low metallicity galaxies). 
All model galaxy parameters were generated randomly, 
and final model galaxy spectra were constructed from the 
BC03 sim ple stellar population (SSP) spectra using the 



Chabrier[ ( J2003) IMF and the Padova 1994 evolutionary 
tracks. 

Each model spectru m was reddened to simulate dust 
in the galaxy using the Chariot & Fall (2000) reddening 
law with reddening E(B — V ) drawn from an empiri- 
cal distribution. This empir ical distribution w as derived 
from Balmer decrement (e.g. Osterbrock & Ferland 2006) 
measurements of star-forming galaxies in the SDSS spec- 
troscopic sample in the redshift range 0.04 < z < 0.10 as 
calculated from the emission line fluxes measured by the 
MPA-JHU team. For galaxies with Balmer decrement 
equal to or below the canonical value a reddening value 



of E{B — V) = was imposed to forbid non-physical 
negative reddening. This distribution effectively acts as 
a prior on the expected amount of reddening in normal 
ga laxies, and is concep tually similar to the one employed 
by [Salimet~aI[ ( [2007[ i. 

We also mspected the use of an effectively flat prior on 
reddening by allowing each model to be fltted for an op- 
timal reddening to match observed photometry values. 
Doing so we found that very blue galaxies with large 
observational error bars were frequently fltted with non- 
physical negative reddening applied to the oldest model 
SFHs. Similarly for red galaxies, each template was fit- 
ted with a large amount of reddening so that the redden- 
ing (and sSFR) estimates were skewed to a distribution 
driven by the difference between the model galaxy color 
distribution and the observed galaxy color. Thus we be- 
lieve this empirical reddening prior is justified in order to 
obtain realistic estimates of reddening in our SN la host 
galaxies. 

Galaxy stellar mass, SFR, and reddening were com- 
puted for each SN la host galaxy as follows. Fluxes for 
our models were synthesized from their SEDs using the 
filter throughputs for GALEX, SDSS, and 2MASS at sev- 
eral steps in redshift {0 < z < 0.2 in steps of Az = 0.005), 
then adjusted according to their assigned reddening. For 
each SN la host, we calculate the fluxes of each model 
SED at the host redshift by interpolating the fluxes com- 
puted at the two nearest redshift steps (we note that our 
choice of flne redshift step resulted in negligible interpo- 
lation error of < 0.005 mag). Each model i was then 
fitted for a simple scaling factor Ai to minimize the x^ 
of its fluxes F\ i matching to the observed photometry 
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We note here that photometric errors (TA,obs for SN la 
host galaxies were padded to ensure physically meaning- 
ful x^ values, as we will justify in detail in Section 3.3 
An example of this fltting procedure is shown in Figure[8 
Each SN la host galaxy property (mass-to-light ratio, 
sSFR, reddening) and its uncertainty was then calculated 
from the probability distribution function (PDF) con- 
structed as the x^ weighted sum of that property across 
all SPS models: 



(p> = E 



-X /2 



(2) 



where the three properties calculated are mass-to-light 
ratio, sSFR, and reddening. In practice we use the x^ 
weights for each model to calculate the weighted cumula- 
tive probability distribution for each quantity, then from 
this calculate median and ±lcr values for the desired 
quantity. 

Galaxy masses were computed using the observed 
g-band absolute magnitude, corrected for the amount 
of dimming Ag corresponding to the fltted reddening 
E{B — V), and the mass-to-light ratio {M^,/L) and its 
error as calculated from the M^,/L PDF. Model speciflc 
SFR (sSFR = the SFR per unit mass) values were calcu- 
lated as the average over the last 0.5 Gyr for each model, 
and the host galaxy fltted sSFR was calculated from the 
sSFR PDF. Similarly, the best flt extinction E{B - V) 
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SNF20070506-006 




Observer Frame Wavelength [A] 



Figure 8. Example SPS photometry fit for the host 

of SNF20070506-006 using three filter subsets; (top) all 
UV/optical/NIR filters, (middle) optical only, and (bottom) only 
the g and i filters. In each panel, the observed photometry points 
are the large green circles, while the purple triangles are the pho- 
tometry points for the single best fitting model SED. The 30 best 
fitting SEDs are shown, color coded so that purple colors indicate 
stronger matches while red colors indicate poorer matches. 



for each host was calculated from the PDF constructed 
from the input model extinction values. 

We explored the possibility of using j-band mass-to- 
light ratios instead of g-hand, since this band ha s for- 
mally smaller scat ter in mass-to-light ratio (see e.g. Gal- 
lazzi fc Bell 2009). We found that for our faintest host 
galaxies z-band yielded a much larger mass error due to 
the blue colors of these galaxies and increased sky bright- 
ness in i-band yielding a larger photometric error than 
g-haiid. For consistency then we chose to use g-hand for 
mass estimates for all our SN la host galaxies. We found 
that the host galaxy mass values derived from i-band 
were consistent with our nominal g-hand mass-to-light 
ratio values, but with larger uncertainties. 

In Table [6] we present the fitted mass, sSFR and host 
reddening E{B — V) fitted for our SNfactory hosts using 
our SPS methods. In Figure [8] we plot a representative 
example fit to photometry for one of our SNfactory host 
galaxies, along with photometry fits using reduced sets 
of filters. The ability of our photometry to constrain 
host galaxy physical parameters was highly dependent 
on the filter sets used, with significant improvement on 



the sSFR and extinction measurements when GALEX 
data was included. We discuss the implications of these 
effects in the following section. 

3.2. The Effect of Uneven Photometric Coverage 

Because our host galaxies do not all have the same set 
of photometric filters, it is vital to ensure that this un- 
even coverage does not bias our results. To this end, we 
tested how well our method was able to recover the input 
model physical parameters when restricted to photomet- 
ric matching with a reduced set of filters. Because a few 
of our SN la hosts were not imaged by GALEX (generally 
due to proximity to bright stars) or were too faint to have 
significant detections in GALEX or 2MASS, we choose 
to examine the effect of having only optical ugriz pho- 
tometry on galaxy parameter recovery. Similarly, since 
our SNIFS observations of faint hosts specifically target 
g and i band photometry, we also examined how well 
galaxy properties could be recovered with only this sin- 
gle color. 

We selected a random sample of 1000 model galaxies 
and perturbed their tabulated magnitudes by values dis- 
tributed normally about zero with a dispersion equal to 
the systematic error values us ed to pad our observational 
photometric data (see Section 3.3). We then applied our 
Bayesian parameter recovery test (excluding each test 
model from the set of models used in the comparison) us- 
ing different subsets of filters: (i) all available filters, (ii) 
optical filters only (no UV/NIR), and (iii) g- and i-band 
only. We show in Figure [9] the comparison of the recov- 
ered mass-to-light ratio, sSFR, and reddening values as 
a function of their input values for each filter subset. A 
representative example of the effects of reduced filter sets 
is also shown for real data in Figure [8J A summary of the 
mean offset (bias) and RMS (dispersion) about the cor- 
rect input values for each parameter for each filter subset 
is presented in Table [7] 

Several important trends are evident in the results of 
this analysis. First, galaxy mass-to-light ratios are re- 
covered with very little bias (< 0.02 dex) regardless of 
the filter set used. This result illustrates the fact that 
galaxy mass-to-light ratios are strongly correlated with 
color, which is set primarily by the mean stellar age of the 
galaxy and only to higher order by the subtleties of the 
galaxy's star-formation history. Furthermore, the color- 
luminosity change due to foreground extinction is only 
slightly dissimilar to the color-luminosity change that ac- 
companies aging in stellar populations. This degeneracy 
is troublesome for accurately estimating the extinction 
and sSFR of a galaxy, but is fortuitous in the recovery 
of mass-to-light ratios from color. Naturally the disper- 
sion in mass-to-light ratio recovery increases as the num- 
ber of filters decrease. Table [7] shows this increase, from 
0.10 dex for the fuU UVOIR filter set, to 0.18 dex for op- 
tical only, to 0.27 dex for gi only. These dispersions are 
satisfyingly small compared to both the range of mass-to- 
light ratios probed by our models (^1.3 dex) and the full 
range of galaxies masses (~5 dex) probed by our SN la 
host galaxy sample. 

The second major result for this analysis is that the 
lack of UV data can severely bias estimates of galaxy 
star-formation intensities (sSFRs) and extinction from 
dust. This effect is particularly acute for galaxies with 
low star-formation intensity (i.e. sSFR < —10.5), where 
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Table 6 

SNfactory SN la Host Galaxy Stellar Population Properties 



SN Name 



log(M./Mo) 



sSFR 
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Figure 9. Model parameter recovery efficiency for a sample of 1000 model galaxies using different subsets of photometric filters. The top 
row shows the recovery of model mass-to-light ratios, specific star formation rates, and reddening by dust when using all filters from the 
UV to NIR. The middle row shows the same results using only optical filters (ugriz), while the bottom row shows the recovery efficiency for 
only two optical filters (g and i). Mass-to-light ratios M/L are in solar units, sSFR values are log(yr~^), and extinction values E{B — V) 
are in magnitudes. Parameter recovery efficiencies are summarized in Table [7] 



the degeneracies between stellar age and extinction be- 
come difficult to disentangle. A galaxy with intermedi- 
ate star-formation intensity and modest reddening (e.g. 
sSFR ~ — 11 and E{B — V) ^ 0.5 mag) has remark- 
ably similar optical colors to an unreddened older stel- 
lar population. Thus galaxies in this sSFR range tend 
to have both their sSFR and extinction levels underes- 
timated. However, galaxies with strong star-formation 



(sSFR > —10.5) generally have successfully recovered sS- 
FRs and extinction, due to the fact that their blue colors 
cannot be mimicked by anything other than young stars. 
As a final remark, we reiterate that these tests were 
done using members of a galaxy model library whose 
parent distribution by construction was the same as the 
model whose properties were being measured. Specifi- 
cally, the reddening for each model was drawn from the 
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Table 7 

Model Parameter Recovery Efficiencies 



Parameter 


Filter Set 


Recov. Eff. 


log(Af./L) 


UVOIR 
Optical 


-0.005 ±0.095 
0.014 ±0.176 
0.024 ±0.265 


log{sSFR) 


UVOIR 
Optical 

g/i 


-0.22 ±0.49 
-0.79 ±0.76 
-1.04 ±1.16 


E{B - V) 


UVOIR 
Optical 

g/i 


-0.012 ±0.016 
-0.065 ±0.053 
-0.087 ±0.085 



same prior distribution it was fitted against, and with the 
same extinction law used as input to the model library. 
Also the simulated measurement errors were distributed 
exactly according to the values used as systematic er- 
rors in our model parameter fits. Thus the results here 
represent the best case scenario in which we fully under- 
stand the underlying distribution of galaxy SFHs, the 
reddening law of the observed galaxies, the distribution 
of extinction values of those galaxies, and the systematic 
uncertainties of our photometric measurements. These 
results then can be considered a systematics floor with 
regard to our ability to recover galaxy physical properties 
from photometry using our SPS method. 

3.3. Validation of SPS Methods 

To ensure the SN la host galaxy properties derived 
with our methods are consistent with those derived by 
other authors, we here examine the results of applying 
numerous SPS methods to the same photometric data as 
a means of validating our method. In order to not bias 
the results of our SN la host galaxy property estimation, 
we conduct our analysis on a training sample of field 
galaxies with photometric data from GALEX, SDSS, and 
2MASS. 

For simplicity (and limitation of computation time) we 
searched the SDSS online database in a small patch of sky 
between -1° < 5 < 1° and 0° < a < 45° (i.e. in Stripe 
82) for galaxies with spectroscopic observations in a red- 
shift range 0.04 < z < 0.10. We then similarly searched 
the GALEX photometric catalog for objects with de- 
tections in both the FUV and NUV filters in the same 
region of sky. Next we searched the 2MASS Extended 
Source Catalog (XSC) in the same sky region. Finally 
we matched the three catalogs together, requiring ob- 
jects be detected with 5" (i.e. the GALEX resolution) 
of the SDSS astrometric position, compiling a total of 
3673 galaxies with reliable detections in all 10 photo- 
metric filters. This set served as our training sample for 
refining our galaxy property estimation method, as well 
as a base sample for comparing the results of employing 
different galaxy SPS methods. 

We note that although these catalog magnitudes do not 
employ matched apertures in the manner of our SN la 
host galaxy photometry (see Section 2.2), they are the 
same magnitudes used in the calculation of galaxy phys- 
ical properties by the MPA-JHU group whose values 
will serve to validate the accuracy of our SPS methods. 
While the advantage of matched photometric apertures 
in galaxy stellar population modeling was not explicitly 



investigated for this work, we found the matched aper- 
ture method especially advantageous for placing upper 
limits on the UV flux in galaxies with low star- format ion. 
The validation sample of SDSS galaxies employed here 
has galaxies spanning the full range of star-formation 
intensity realized in our SPS models, so the choice to 
use catalog photometry has not excised low-SFR galax- 
ies from our validation sample. 

Because our galaxy property (and uncertainty) esti- 
mation relies on appropriate weighting of the models in 
our library, it is vital that our photometric uncertainty 
estimates be accurate. We inspected this by examining 
the distribution of best x^ values for our training sam- 
ple of SDSS galaxies and calculating its agreement with 
the expected distribution of this quantity. We find that 
use of the nominal catalog photometric errors results in 
extremely large x^ values. We thus add extra systematic 
errors in quadrature with the formal photometric errors 
of our measurements as follows: 0.052 and 0.026 mag for 
GALEX F UV and NUV filters, respectively (c.f. [Salm 
et al.|2007[ ); 0.05 0.02, 0.02, 0.02, and 0.0 3 mag for^DSS^ 
ugriz tilters (c.f. Blanton & Roweis 2007); and 0.05 mag 
for each of the 2MASS and UKIDSS JHK filters. The 
distribution of best x'^ values using these padded pho- 
tometric error bars is in better agreement with the ex- 
pected x^ distribution, so we employ them in the galaxy 
property estimation for our SN la host galaxies. 

We note that the final y^ values are still slightly higher 
than predicted by the x distribution, indicating some 
residual disagreement between our models and the data. 
This is likely due to the fact that the true underlying 
distributions of galaxy SFHs, reddening laws, and extinc- 
tion values are not fully known quantities. Similarly, the 
best stellar evolutionary tracks, IMF(s), and even stellar 
SEDs may still require refinement. While these quan- 
tities are of interest for proper galaxy modeling, their 
study is the subject of a large and rich field of inquiry 
and thus must be beyond the scope of this work. Our 
choices of SFH, reddening, and metallicity priors as well 
as IMF, stellar evolutionary track, and systematic photo- 
metric uncertainties, are all based on reasonable assump- 
tions which have even been employed in previous studies. 
Ultimately our models show acceptable agreement with 
real data and thus we consider them to be good models 
for comparing to data. 

To confirm that our Bayesian method is recovering the 
correct galaxy physical properties, we compare the val- 
ues we derived for our training sample of SDSS galaxies 
to the values derived by previous authors. To this end 
we use the stellar mass and SFR values calculated by the 
MPA-JHU group for these same galaxies. In the left hand 
panels of Figure [lO] we show the values we derive com- 
pared to their values. We find our stellar masses to have 
excellent agreement, with a mean and RMS difference of 
Alog(A/) — 0.063 ±0.105 and typical galaxy mass errors 
of '^\og{Ai} — 0.16. Their values use the same SPS mod- 
els (BC03) as us with the same IMF ( |Chabrier|2003| , as 
well as model galaxy SFH parameters drawn from a sim- 
ilarly constructed prior and a similar Bayesian weighting 
scheme, so such tight agreement is to be expected. 

Our SFR values show a bit more dispersion compared 
to the MPA-JHU values, with Alog(5Fi?) = -0.4 ± 1.2. 
This is not unexpected for several reasons. The first and 
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most important is that the MPA-JHU SFR values are 
based on a combination of the SFR calculated from Ha 
flux within the SDSS fiber and color-based SFR mea- 
surement for light outside the fiber, which is trained on 
the Ha fluxes and colors of fiber spectra. The second 
cause for SFR discrepancies is the fact that the Ha- 
trained SFR values for the SDSS sample are sensitive 
to the most recent star-formation and may not correlate 
exactly with our UV-driven SFR values averaged over 
the last 500 Myr of the galaxy SFH. These effects will be 
discussed in detail in Section [5l 

Next we compute mass and SFR estimates for our 
training sample of galaxies using the code ZPEG (Le 
Borgne & Rocca-Volmerange 200_2). This code was em- 
ployed by several recent authors ( Kelly et al.|2010 Sulli- 
van et aL]|2010 Lampeitl et al.pUlOp to examine trends 
in stretch- and coior-corrected SiN la Hubble residuals 
with the masses of their host galaxies, and thus we want 
to ensure our values are on a consistent scale with theirs. 
We used our catalog photometry for our training sam- 
ple of galaxies and computed masses and SFRs using 
ZPEG with the standard set of templates included with 
the ZPEG code. The comparison of the values derived 
by the MPA-JHU groups and those derived with ZPEG is 
presented in the right panels of Figure [TOJ which shows 
that the stellar mass values derived with ZPEG are highly 
consis tent with those de r ived b y the independent analy- 
sis of Kauffmann et al. (2003). The comparison of the 
stellar mass values frorn ZPEG with our own values also 
showed similar agreement. We note however that the dis- 
crete set of input galaxy evolutionary scenarios in ZPEG 
results in nearly discrete values of fitted sSFR. 

We found that the formal ZPEG mass and sSFR er- 
ror bars are significantly smaller than the values derived 
with our method. This is because ZPEG calculates errors 
on these quantities based on how well the input pho- 
tometry matches the single evolutionary scenario that 
fits best, but has no prescription for how uncertain the 
galaxy SFH might be. Our method implicitly incorpo- 
rates this estimation by performing a weighted sum over 
a wide range of galaxy SFHs, so we believe that our error 
bars more accurately reflect the systematic errors in the 
galax y mass and SFR estimation (see also Taylor et al. 
|2011[ for a thorough discussion on the superior results 
achieved with Bayesian galaxy parameter measurements 
compared to using the "best-fit" template from a discrete 
set of models) . 

Thus we conclude that our mass values are appropri- 
ate for use in SN la cosmological analyses and a re con- 
sistent with other methods including ZPEG and Kauff- 
mann et al. (2003). Indeed, since the central values of 
all these metnods are so consistent, it is highly unlikely 
that trends of SN la Hubble residuals with host galaxy 
mass can be measured to higher precision by employing 
more sophisticated SPS methods. However, if these host 
galaxy mass values should ultimately be used as SN la 
luminosity correction parameters, then the estimation of 
their error bars must be accurate in order to be properly 
accounted for in the SN la cosmology error budget. Thus 
we propose that our Bayesian method is more realistic as 
it accounts for galaxy modeling systematics where some 
other methods do not. 
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Figure 10. Left Panels: Comparison of galaxy stellar masses 
(top) and specific star-formation rates (bottom) derived using the 
Bayesian SPS method described here versus the values computed 
by the MPA-JHU SDSS team for a sample of galaxies in SDSS. 
Right Panels: Same as left, but comparing the ZPEG values to those 
of the MPA-JHU group. 



3.4. The Importance of GALEX UV Photometry 

We further investigated the importance of including 
UV photometry from GALEX by repeating our galaxy 
photometry fits without UV data, both for the training 
sample of SDSS galaxies and our sample of SN la host 
galaxies. In particular we focused on the change in the 
estimated galaxy sSFR when UV data was not included 
in the photometry fits. 

For the sample of ~ 3700 SDSS galaxies in our train- 
ing sample, we found that the sSFR values calculated 
from fitting only optical data were on average 0.22 dex 
lower than the values obtained when utilizing GALEX 
data. This is in good agreement with Figure [9| where we 
found that the lack of UV data resulted in many inter- 
mediate sSFR galaxies with modest dust being mistaken 
for low sSFR galaxies with no dust. Furthermore, the 
mean uncertainty in the estimation of sSFRs went up 
from 0.2 dex to 0.72 dex when UV data was not included 
in the SPS fit. Thus we see that UV data is critical 
for both obtaining the correct value of a galaxy's star 
formation intensity as well as reducing its measurement 
uncertainty. 

We also specifically investigated the effect of neglect- 
ing UV data when fitting the photometry of massive el- 
liptical (passive) galaxies. To investigate this portion of 
parameter space, we isolated the subset of ^ 950 galax- 
ies in our training sample with log(M*/M0) > 10.5 and 
\og{sSFR) < —11.5. For these galaxies we found the 
sSFR values fitted from only optical photometry were 
0.28 dex higher than their values obtained from fitting 
the full UVOIR data. These galaxies correspond to the 
models in the lower left of the middle panels of Figure [9j 
where the sSFR is over-estimated in extremely low sSFR 
systems. These systems behave in the opposite fashion to 
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their intermediate-sSFR counterparts, as their red colors 
are due to extremely old stars but the prior on dust red- 
dening forces both their fitted sSFR and dust reddening 
values to be too high. UV data breaks this degeneracy 
and allows the extremely low sSFR and low dust extinc- 
tion to be correctly recovered. 

These trends observed in the photometry fitting of 
both the models and SDSS galaxies were also observed 
(albeit with smaller statistics) in the SNfactory SN la 
host galaxy sample. For the 394 hosts with GALEX 
photometry, the mean sSFR uncertainty increased from 
0.27 dex to 0.43 dex when the UV data was not in- 
cluded. We found the GALEX UV data to be critical 
for correctly determining the sSFR of our massive el- 
liptical hosts. For the 66 SNfactory host galaxies with 
log(M*/Af0) > 10.0 and log{sSFR) < -11.0 (note we 
softened the cuts to increase the sample size), the sSFR 
values were on average 0.40 dex higher when GALEX 
UV data was excluded in the photometry fit. Of those 
hosts, 35 have GALEX data with lower than 5a detec- 
tions, illustrating that even modest detections or upper 
limits of flux in the UV provide critical information for 
the determination of SN la host galaxy star formation 
intensity. 

4. SN la HOST GALAXIES AND THE GALAXY 
MASS-METALLICITY RELATION 

The level of agreement of SN la host galaxies with the 
normal galaxy mass-metallicity (MZ) relation can pro- 
vide important insight into preferred SN la progenitor 
environments. Discrepancies between the SN la MZ dis- 
tribution and that of normal galaxies could potentially 
indicate metallicity preferences for SNe la, which would 
have important implications for high redshift SN la sur- 
veys. Alternatively, disagreement with the MZ relation 
could have other interpretations, as was the case with 
long-duration gamma ray burst host galaxies. 

Some recent studies of the host galaxies of long- 
duration gamma ray bursts (LGRBs) found that they 
tended to have systematically lower metallicit ies than 
those predicted by fiducial galaxy MZ relation (Modjaz 
et al.|[2008 Levesque et al.pOlOJ. Initial interpretations 



of this trend speculated on a preference for lower metal 
licity environments in the production of LGRBs. The 
key insight, however, came from considering the effect 
of galax y star-formation rate ( SFR) on the galaxy MZ 
relation (Mannucci et al. 2010). Accounting for this ef- 
fect, it was found that LGRB hosts indeed agreed with 
the SFR-adjusted MZ relation (or equivalently the M- 
Z-SFR relation) but merely appeared in the region of 
galaxy parameter spa ce populated by the most intensely 
sta r-forming galaxies ([Kocevski fc West]|2011 Mannucci 



|et al.||201l] ). Thus this trend showed the preference tor 
LGRBs to form in very young stellar environments. 

The SN la host galaxy agreement with the MZ rela- 
tion has been an implicit assumption of previous authors 
who interpreted SN la brightness trends with host galaxy 
stellar mass in terms of SN la progenitor metallicity. The 
SNfactory sample is ideal for testing this assumption, as 
our untargeted search found SNe la in an unbiased sam- 
ple of host galaxies spanning a large range (nearly 5 dex) 
of stellar masses. In this Section we present our method 
for inspecting the consistency of SN la host galaxies with 
the galaxy MZ relation and the results from the hosts of 



SNe la discovered by SNfactory. 

4.1. SNfactory SN la Hosts and the MZ Relation 

The correlation of galaxy luminosity and stellar mass 
with metallicity ha s been known for several decades 
(Lequeux et al. 1979 ), but has been quantitatively refined 
only recently with the adven t of major galaxy spectro- 
scopic surve ys at low (SPSS York et al. 2000) and in- 
termediate (Zahid et al. 2011[ ) redshitts. Of particular 
interest for this work is the correlation of galaxy stel- 
lar mass with gas-phase metallicity, which for simplicity 
we will refer to simply as "metallicity" in this Section. 
For SDSS the MZ relation was studie d by the MPA-JHU 
SDSS team in |Tremonti et al.|p004l T04) for the fourth 
SDSS data release and subsequently for future data re- 
leases. T04 found that for a sample of «45,000 galaxies, 
gas-phase metallicities followed a tight relation in the 
stellar mass range of 8.5 < log(Af*/Af0) < 11.0 with a 
dispersion of about 0.1 dex at high stellar masses. The 
dispersion in the MZ relation increases at lower stellar 
mass, up to about 0.3 dex at log(M*/M0) = 8.5. 

For this analysis we wish to inspect how much SN la 
hosts deviate from the fiducial MZ relation and whether 
those deviations are consistent with the observed disper- 
sion in the MZ relation. To do so we use derived stellar 
masses and metallicities fro m the MPA-JHU SPSS team 
analysis of the SDSS DR7 ( [Abazajian et al.||2009[ ) data. 
They derive galaxy stellar masses from broadband pho- 
tometry using the stellar population synthesis library of 
Kauffmann et al. ( 2003) , and calculate gas phase metal- 
licities from emission line fluxes according to the method 
outlined in T04. To facilitate the appropriate compari- 
son, we use the SNfactory host stell ar m asses and metal- 
licities derived in Sections [2?2] and [2!3l Th e MPA-JHU 
masses are computed with a Kroupa ( 20011) IMF , which 
has a negligible mean offset from tne Ghabrier (2003) 
IMF employed in our methods, though there exists some 
small scatter. As stated above, our metallicity values 
have all been converted to the T04 scale, so our masses 
and metallicities can be compared directly to those de- 
rived for SDSS galaxies by the MPA-JHU group. The 



full SNfactory host MZ diagram is shown in Figure 11 

In order to assess the agreement of SNfactory host 
masses and metallicities with the SDSS MZ relation, 
we first compare the observed metallicity of each SN la 
host galaxy with the value predicted by the MZ rela- 
tion for its observed mass. In practice, we calculate 
the weighted mean of the metallicities of all neighbor- 
ing (in mass) SDSS galaxies with masses M^, weighted 
by their distance from the observed host mass Mhost (i-e. 
exp[-x- /2] where xf = {{M, - Mhost)/crA/r)^ and cta/ is 
the uncertainty in the SN la host mass) with proper ac- 
counting for the number of SDSS hosts as a function of 
mass. Thus for each SN la host we can calculate the 
difference between its observed metallicity and that pre- 
dicted from the MZ relation as AZ = Zhost — ZmZi with 
an uncertainty equal to the quadrature sum of the host 
metallicity measurement error and the dispersion of the 
MZ relation at that host mass (i.e. the RMS of the metal- 
licity values of its stellar mass neighbors from SDSS). 

Performing this calculation for all 130 SN la hosts 
in the SNfactory sample in the stellar mass range over 
which the MZ relation is well populated by SDSS, i.e. 
the aforementioned 8.5 < log(Af*/Af0) < 11.0, we find 
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Figure 11. Left: Location of SNfactory iiost galaxies in tlie MZ plane using final T04 metallicities. The grey background is a density 
plot of the galaxies in the SDSS DR7 sample analysis from the MPA-JHU team. Top Right: PP04 N2 metallicities for all SNfactory host 
galaxies. Those hosts whose final metallicities were combined using this method are shown as solid green circle, while those with final 
metallicities calculated from iJ23 are shown as open red squares. Bottom Right: Same as top right, but showing R23, metallicities for 
all SNfactory hosts, with filled red squares showing those hosts whose finally metallicities came from iJ23 and open green squares those 
calculated from PP04 N2. 



the weighted mean (and RMS) deviation of SN la host 
metalhcities from the MZ relation to be: 



(AZ) =0.0007 ±0.1522 



(3) 



For this number of hosts, the error on the mean is 
0.012 dex. If we rephrase the MZ deviation in terms 
of puh values we find a mean and RMS deviation of: 



AZ\ 



= 0.00 ±1.12 



o-z 



(4) 



Thus we can see that SN la host galaxy metallicities are, 
on average, very consistent with the values predicted for 
their host masses by the galaxy mass-metallicity relation. 
The fact that the RMS of our pull values is close to 1 
implies agreement of SN la hosts with the MZ relation 
for not only the average metallicity values, but also the 
observed dispersion. This result also implies that our 
measurement error bars are not misestimated. 

In the right panels of Figure [TT] we plot the SNfactory 
host galaxy MZ relation for the i?23 and N2 metallic- 
ity metrics from whi ch the final T04 metallicities were 
derived (see Section 2.3.3). These plots illustrate the 
aforementioned sensitivity failure of the Nil line at low 
metallicities, as well as the double-valued nature of the 
i?23 metric. As can been seen from these plots, nearly 
all of the SN la host metallicities used in the MZ agree- 
ments tests of this Section come from the N2 metric, 
and show a similar agreement to the MZ relations of the 



SDSS sample using that metric. 

To calculate each individual host galaxy's deviation 
from the SDSS MZ relation in further detail, we derive 
the metallicity cumulative distribution function (CDF) 
at each value of stellar mass, again using the weighted 
metallicities of each host's neighbors in stellar mass. The 
host is then assigned a score corresponding to where its 
metallicity is placed in the CDF of metallicities at its 
mass. This principle is illustrated in the left panels of 
Figure [T2j 

Thus for each SN la host galaxy, we have a measure of 
where its metallicity lies in the distribution of metallici- 
ties at its stellar mass, which we will call its MZ agree- 
ment score. If SN la hosts perfectly obey the MZ relation, 
then the ensemble distribution of these scores should be 
distributed uniformly bet ween and 1. We show in the 
right panels of Figure [12] this distribution of MZ agree- 
ment scores for the 130 SNfactory SN la host galaxies 
whose mass falls within the aforementioned range. From 
this histogram we can see that the scores are relatively 
uniform. In the right panel of the same Figure, we plot 
the cumulative distribution function of the MZ agree- 
ment scores as compared to a line of unit slope (i.e. 
the CDF for a flat distribution). We can see from this 
plot that indeed our distribution is very close to a flat 
(uniform) distribution, and the cumulative distribution is 
reasonably close to unity. This would imply that not only 
are the mean and RMS metallicity deviation for SN la 
hosts consistent with the MZ relation, but the shape of 
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Figure 12. Top Left: Example of method for calculating the MZ agreement score for each SN la host in the SNfactory sample. The blue 
and red bars correspond to the itlcr mass values for two hosts, with the white circles showing their mass and metallicity values. Middle Left: 
The (unweighted) histogram of metallicities within ±0.05 dex of each host mass. Bottom Left: The (weighted) cumulative distribution 
function (CDF) for metallicities at each host's mass. The score for each host is the intersection of its metallicity value (vertical lines) with 
the metallicity CDF at its mass. Top Right: Distribution of MZ scores for SN la hosts from SNfactory. Bottom Right: CDF comparison 
of the SN la host MZ score distribution to a uniform (flat) distribution. 



their distribution is also similar. 

We note that 30 hosts in our sample have host masses 
below the lower mass limit over which the MZ relation 
is well defined. We measure the RMS of the metallic- 
ity values of SNfactory hosts at the low mass end to be 
0.32 dex, which is consistent with the 0.3 dex dispersion 
me asured by T04 and the 0.4 dex dispersion measured 



by Zahid et al. (2012) at comparably low masses. In- 
terestingly there appears to be a paucity of hosts below 
12 + log{0/H) = 7.7, near the low me t allicity thresh- 
old predicted by Kobayashi & Nomoto (20091. While 
interesting, we caution that proper interpretation of this 
result requires a complete accounting of both observa- 
tional incompleteness and SN search efficiency. A careful 
treatment of this is underway and will be presented in a 
future SNfactory paper. 

5. STAR-FORMATION ACTIVITY OF SN la HOST 
GALAXIES 

In a similar spirit to the analysis of Section [4j we wish 
to investigate whether the star- format ion rates of SN la 
host galaxies are distributed similarly to the values one 



would measure from a comparable sample of normal field 
galaxies. Galaxy SFRs are not tightly correlated with 
stellar mass as was true for gas-phase metallicities, so 
here we will examine whether the distribution of star- 
formation indicators for the SNfactory SN la host galaxy 
sample is similar to the same distribution for a compa- 
rable sample of field galaxies. 

This analysis begins with the expectation of finding 
subtle differences between SN la hosts and field galax- 
ies driven by the SN la delay time distribution (DTD) 
- the distribution of times from formation of the pro- 
genitor system to explosion as a SN la. SN la rates 
show a d ependence on b oth host gal axy SFR and stel 
lar mass (Mannucci et a l. 2005; Scan napieco fc Bildsten 
f^OOF; Mannucci et al. 2006; S ufiivan et aL|20Ut)[|Aubourg 
let al 2008; Smith et al 2012"| , and studies of the SN la 
D'l'D indicate a larger traction of young (< 1 Gyr) pro- 



genito rs CTota ni et al.||2008l [Brandt et al. [20101 [Barbary 
|et al.| [2012; M aoz et a!l .||2012p. I'herefore, we begin this 



comparison with the expectation that our SN la host 
galaxies may show enhanced star formation. 
In order to compare SN la host galaxy SFRs to a sam- 
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pie of typical galaxies, we require an appropriate field 
galaxy sample - ideally one with stellar masses and SFRs 
measured from photometry in a manner comparable to 
our own methods, as well as spectroscopic measurements 
of galaxy Ha emission. While there are many larg e well- 
studied sam ples that come close (e.g.,_ NFGS (.Jansen 



photome tric SFRs and S FRs calculated from Ha flux 
( |1998[ ) relation (Figure [?]). 



with the Kennicutt 



We now turn to t 



et al.|2000| ; IIHUGS ( [Kennicutt et al.||2008[ [Lee et al. 
2Ulip ), the galaxy catalog that currently provides the 
closest match to our needs is again the MPA-JHU value- 
added SDSS catalog. We will use this as our baseline for 
establishing the SF properties of a field galaxy sample, 
but in what follows we will describe a number of system- 
atic limitations that do not yet allow us to exploit the 
full statistical power of our SN la host galaxy sample. 
Thus, the aim of this Section will be to test whether the 
star-formation activity of SN la host galaxies differs from 
that of a comparable sample of field galaxies at a level 
significantly above the systematic uncertainties. 

Accordingly, we proceed by first discussing some key 
differences between the MPA-JHU method and our own 
tech niqu e for estimating SFRs from photometry in Sec- 
tion |5.1[ Then in Section [5.2| wc discuss the importance of 
appropriately weighting the SDSS spectroscopic sample 
to account for SFR-dependent Malmquist bias. Finally 
in Section |5.3| we present the results of our SN la host 
galaxy comparison to SDSS field galaxies. 

5.1. Differences in the SDSS and SN factory SFR 
Measurement Methods 

Photometric SFI ls for the MPA-JH U SDSS catalog 
were presented in Brinchmann et aT] ([2004, hereafter 
B04). These SFRs comprise two components: the first 
calculated directly from the spectrum obtained from the 
3" fiber placed on the center of the galaxy, and the sec- 
ond calculated from the flux and colors of the galaxy 
light outside the fiber. 

The fiber SFR values in B04 were calculated from a 
Bayesian technique similar to the photometric mass es- 
timates applied to the measured emission line fluxes in 
the galaxy spectrum. These measured line fluxe s were 
compared to a suite of emission line models from |Char-| 
lot & Longhetti (2001) which spanned various values of 
metallicity, ionization parameter, and reddening by dust 
(note these are the same models and comparison meth- 
ods used to calculate metallicities in T04). B04 found 
the flnal SFRs to be driven primarily by the flux in Ha, 
but with a scaling parameter that varied subtly from the 
Kennicutt ( 1998 1 value depending on the best flt metal- 
licity and ionization parameter. 

SFR values for the galaxy flux outside the SDSS fiber 
were calculated in B04 using a color-color SFR grid 
trained on the spectroscopic data. Specifically, B04 cal- 
culated the SFR per unit luminosity in i-band within the 
SDSS fiber for all galaxies in their sample, then measured 
the mean SFR/ Li in small bins in g — r vs. r — i color 
space. For each galaxy, the g — r and r — i colors of the 
light outside the SDSS fiber was measured and the mean 
SFR/ Li for the corresponding color-color bin applied to 
the j-band fiux outside the fiber. 

Thus the B04 SFR values for the SDSS galaxy sample 
were trained to show tight agreement with SFRs mea- 
sured from Ha flux (indeed they are partly composed 
of the Ha flux-based SFR from within the spectroscopic 
fiber). This explains the tight agreement between their 



le systematic differences oetween the 
B04 SFR technique and that employed by us (Section [3| 
for SN la host galaxies. Our technique measures tne 
SFR of a galaxy from its average sSFR over the past 
500 Myr, and is strongly driven by UV flux. This in- 
herently probes a different timescale in the galaxy SFH 
than the Ha-driven SFRs of B04, which probe the very 
recent (^10 Myr) part of the galaxy SFH. There is cause 
to believe that these SFRs should agree on average but 
have some scatter due to the variety of SFHs reali zed in 
nature . This was borne out by the work of ,Salim et al.| 
( 2007 1 - who performed a SED matching technique sim- 
ilar to, and the inspiration for, our own technique - who 
found that the photometrically measured SFRs agreed 
with the Ha-driven SFRs with a scatter of 0.5 dex. Their 
SFRs were calculated over the last 100 Myr in their mod- 
els, and the scatter is likely a product of the diversity of 
galaxy SFHs. Similar scatter (perhaps even larger) is 
to be expected when comparing our 500 Myr averaged 
SFRs to the Ha-driven ones of B04. 

Differences between our SFRs and those derived with 
the B04 technique may also arise due to the fact that 
their values are derived (in part) from optical photome- 
try. This could yield a sS FR-correlated bias as found in 
our models in Section [3.2[ and may explain t he st ructure 
in the SFR comparison presented in Figure |10[ These 
concerns, along with the SFH diversity scatter described 
above, limit the power of direct comparison between SFR 
values derived with these differing method. Our objec- 
tive in this Section then will be to look for strong devi- 
ations in the SFR behavior of SN la hosts which cannot 
be explained by these systematic differences, while a de- 
tailed quantitative comparison will be best suited for a 
future analysis where consistent methods can be applied 
to SN la hosts and a large held galaxy sample. 

5.2. Weighting SDSS Galaxies for Selection Bias 

The objective of this section is to the compare the dis- 
tribution of SN la host galaxy SFRs, as traced by several 
indicators, to a comparable field galaxy sample from the 
MPA-JHU SDSS catalog. Unlike our previous analysis 
of the galaxy MZ relation, this SFR comparison must ac- 
count for observational incompleteness resulting from the 
magnitude limit (m^ — 17.77) of the SDSS spectroscopy 
survey. More intensely star-forming galaxies will be bluer 
and more luminous (i.e. have a lower mass-to-light ra- 
tio) than less intensely star-forming galaxies of the same 
mass. This means that the raw demographics of galaxies 
in the SDSS spectroscopic survey will be biased toward 
more strongly star-forming galaxies. This effect was not 
important in our consideration of the galaxy MZ relation 
(Section l4|) since metallicity does not strongly affect the 
galaxy luminosity, but it is of critical importance for the 
analysis of this Section. 

To mitigate the effect of luminosity bias, we calculate 
for each galaxy in the MPA-JHU sample the maximum 
possible redshift zum at which it would have been in- 
cluded in the SDSS spectroscopic sample, then calcu- 
late the effective volume of the SDSS sample in which 
these galaxies would be found {V ex (zfj^ — z'^in) where 
Zmin = 0.02 is the low-redshift cut for our sample - see 
below). Each galaxy then is assigned a weight accord- 
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ing to its effective fraction of the total survey volume 

iylVrnax wherC Vmax OC (z^^^ - zi^,„)). 

This Y lYmax weighting scheme is complicated by the 
fact that the SDSS spectroscopic sample includes some 
galaxies whose luminosity is below the nominal magni- 
tude limit of the survey. These are typically low mass 
galaxies whose blue colors and small physical sizes made 
them prime candidates to be high-redshift quasars, which 
were spectroscopically targeted below the nominal mag- 
nitude cut. Indeed, below galaxy mass of log(M*/M0) ~ 
8.5 there are almost no galaxies which are more lumi- 
nous than the nominal SDSS spectroscopic magnitude 
limit given our lower redshift cut. For completeness we 
wish to include these galaxies in order to characterize 
the SFR activity at low galaxy mass scales, but we must 
weight them differently than galaxies brighter than the 
nominal SDSS magnitude limit. Upon inspection of the 
magnitudes of these galaxies, we found that most of these 
galaxies obeyed an apparent upper magnitude limit in r- 
band of m^ = 20. We therefore weighted all galaxies 
above the nominal SDSS spectroscopic magnitude limit 
of TOr = 17.77 using a V iVmax weight calculated using 
m^ = 20 as the survey magnitude limit. 

Since galaxy star-formatio n intensity (name ly sSFR), 
correlates with galaxy mass (Salim et al. 2007), we con- 
struct the distributions of star- format ion indicators from 
the SDSS sample using final weights that yield a weighted 
galaxy mass distribution equivalent to that of our SN la 
host galaxy sample. First, we employed redshift cuts 
of z^in = 0.02, to mitigate the effects of low fiber fill- 
ing fractions for massive galaxies at low redshift, and 
Zmax = 0.10, above which the SDSS spectroscopic sur- 
vey is heavily biased towards massive luminous galax- 
ies. V/Vmax volumetric weights were calculated for each 
galaxy using these values and the observed galaxy lu- 
minosity in r-band as described above. We then con- 
structed the volume-weighted galaxy mass distribution 
for this sample and calculated a mass-based weighting 
function to force the final weighted mass distribution to 
match that of the SN la h ost s ample. 

In the analysis of Section [5.3[ we will examine SF indi- 
cators for two subsets of the SNfactory sample: (i) those 
with SFRs (and sSFRs) measured from photometry, and 
(ii) those with Ha SFRs measured from spectroscopy. 
Since these are two distinct samples with different mass 
distributions, we calculated two mass-based weighting 
functions for the SDSS sample to make the appropriate 
comparisons to the respective samples. 

5.3. Comparison of SN la Host Galaxy SFRs with 
SDSS Field Galaxies 



Figure 13 presents a graphical comparison of star- 
formation indicators for the SNfactory SN la host galaxy 
sample compared to the appropriately weighted distribu- 
tion of those same quantities from the MPA-JHU SDSS 
galaxy sample. The top left and right panels show sSFR 
and SFR measured from photometry, while the bottom 
left and right panels show the equivalent width of the 
Ha line and total galaxy flux in the Ha line. For this 
Figure the sSFR and SFR values from the MPA-JHU 
catalog have been randomly perturbed by values gaus- 
sianly distributed about with a dispersion of 0.5 dex, 
in order to simul ate the effect of the UV-Ha SFR disper- 
sion observed by Salim et al. (2007) which we suggested 
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Figure 13. Top Panels: Photometric estimates of sSFR (left) 
and SFR (right) for SNfactory SN la host galaxies (solid blue his- 
tograms) compared to the weighted distribution of those quanti- 
ties from a comparable sample of SDSS galaxies (dashed magenta 
curves). Bottom: Spectroscopic measurements of the equivalent 
width of (left) and total flux in (right) the Ha line, with curves the 
same as top panels. 



Table 8 

Properties of SF Indicator Distributions 



Quantity 



Sample 



-la Median -l-la" 



Value 



Value Value 



log(sSFR) 


- phot 


SDSS 


-10.61 


-9.95 


-9.29 






SNf 


-10.89 


-9.93 


-9.18 


log(SFR) - 


phot 


SDSS 


-0.86 


0.00 


0.82 






SNf 


-1.12 


-0.18 


0.59 


EW(Ha) 




SDSS 


13.6 


24.2 


47.0 






SNf 


13.9 


30.0 


46.9 


log{F{Ha)) 


SDSS 


40.21 


40.86 


41.42 






SNf 


39.92 


40.72 


41.27 



above (Section 5.1 ) is likely to arise from SFH diversity. 
We note that the sSFRs for our SN la hosts are averaged 
over the last 500 Myr of the model SFHs and thus by con- 
struction have a maximum value of log{sSFR) = —8.7 
corresponding to the entirety of the galaxy mass being 
formed in this time interval. 

In Table [8] we summarize the properties of the star- 
formation indicator distributions for the SDSS galaxy 
sample and the SNfactory SN la host galaxies (we note 
that for this comparison we restrict the sample only to 
SNe la discovered by SNfactory). We present the me- 
dian and ±1ct values for the two samples, but due to the 
systematic differences i n ph otometric SFR measurement 
techniques (see Section 5.1) and t he uncertain impact of 
aperture effects (see Section 2.3), we purposfuUy avoid 



power statistical tests such as a Kolmogorov-Smirnov 
test because they cannot effectively account for these 
subtleties. 
From this analysis we can make some general state- 
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ments about the consistency of SN la host galaxy star- 
formation activity compared to a comparable sample of 
field galaxies. The differences between photometric SFR 
and sSFR distributions is much smaller than the system- 
atic uncertainties in their measurements, so we detect 
no significantly unusual star- format ion activity in SN la 
host galaxies as measured from photometry. 

Interestingly, the EW{Ha) distribution for SN la host 
galaxies appears to have a higher median value than the 
SDSS field galaxy sample. This could perhaps indicate 
that SNe la have a slight preference for more strongly 
star-forming environments. This would be consistent 
with recent studies of the SN la delay time distri bution 
et al.|2011[|Barbary et al.|2012[ ) -5 
d at young stellar ages. Quantit 



(DTD - e.g. 
find a DTD 



Maoz I 



which 



peaked at young stellar ages. Quantitative 
results from this analysis are limited by the uncertain im- 
pact of aperture effects, but EW{Ha) provides the most 
interesting distinction between the SN la host galaxy and 
normal field galaxy samples. 

In summary, SN la host galaxies show mostly normal 
star-formation activity given the current limitations on 
measurement of SFRs in our host galaxy sample and the 
field galaxy sample from SDSS. A more precise quantita- 
tive comparison will likely require a field galaxy sample 
whose properties are derived in a consistent way as the 
SN la hosts, and whose selection function is very well 
understood. Since SN la hosts have a higher relative 
fraction of low mass galaxies compared to the magnitude- 
limited SDSS sample, this comparison would also benefit 
from a field galaxy sample which is better sampled at 
the low mass end. For now we find modest agreement 
between the star-formation activity of SN la host galax- 
ies and a comparable (in mass) sample of field galaxies. 

6. SUMMARY 

Here we revisit the major contributions of this work. 
Specifically, we summarize the properties of the SNfac- 
tory SN la host galaxy sample compared to other SN la 
host samples, and illustrate that the SNfactory SN la 
sample found SNe la in environments missed by other 
surveys. We then discuss the advantages of our galaxy 
stellar population synthesis (SPS) technique, as well as 
its possible limitations. Finally we discuss the impli- 
cations of our findings that SN la host galaxies exhibit 
similar metallicities and star-formation behavior as com- 
parable samples of normal field galaxies. 

6.1. Advantages of the SNfactory SN la Sample 

The SNfactory SN la sample whose host galax- 
ies are presented here represents a sample of 400 
low redshift {z < 0.1) SNe la from an untargeted 
search, which we believe is the best available repre- 
sentation of the full population of SNe la and their 
host galaxies in the local universe. A representa- 
tive statistic is the range of stellar masses probed by 
the sample of SNe la discovered by SNfactory, which 
spans from log(M*/Af0) = 6.9 (SNF20080910-007) to 
log(M*/M0) = 11.7 (SNF20050926-002), with nearly 
half (46%) of our host galaxies having stellar masses be- 
low 10^^ Mq. As a comp arison, the compil ation of local 



We also used optical spectra to measure the gas phase 
metallicities for a large fraction of the SNfactory host 
galaxy sample, and showed that our hosts span nearly 
2 dex in metallicity. This is especially advantageous for 
the construction of SN la spectral templates for applica- 
tion at high redshift where metallicities and stellar ages 
will be on average lower than in the low redshift universe. 
Indeed, a measuremen t of the galaxy ma ss-metallicity 
relation at z ~ 0.8 by |Zahid et al. (20111 showed that 
the MZ relation has shitted downward by approximately 
0.15 dex, less than a tenth of the full metallicity range 
probed by our sample. Thus we can state that our ability 
to construct appropriate SN la templates at low redshift 
will not be hindered by cosmic chemical evolution since 
the range of metallicities probed locally far exceeds its 
average change at high redshift. 

6.2. SN la Host Galaxies and SPS Techniques 

In Section [3] we presented our Bayesian technique for 
deriving SN la host galaxy masses and specific star for- 
mation rates from multi-band photometry. This method 
of applying stellar population synthesis (SPS) techniques 
to SN la host galaxies differs somewhat fr om those used 
by other authors studying SN la hosts (e.g . Sullivan et al 



2006 



2010 



Neill et al.||2009[ jKelly et al.|[20To| |Sulliyan eTal 
Lampeitl et al.l 2010p. We made several key ob- 



SN la host galaxies from Neill et al. ( 2009 1 spanned the 
mass range 8.6 < \og{M^/MQ) < UM^ith only 15% 
of their hosts having stellar masses below lO^^M©. 



servations about our technique which are important for 
future comparisons of SN la properties to the properties 
of their host galaxies. 

First, we noted that SPS fits to galaxy photometry re- 
covered the star-formation rate well only if UV data are 
available. Because of the degeneracy between redder col- 
ors in old stars and reddening by foreground dust, optical 
data alone cannot effectively distinguish between these 
two effects. Instead, both the galaxy sSFR and dust are 
underestimated at intermediate galaxy ages with modest 
reddening by dust, where the SPS recovery method mis- 
takes the dust for older stellar ages. Thus our method 
is optimum when rest-frame UV photometry is available. 
For high-redshift ground-based surveys such as SNLS and 
DES, observer- frame optical photometry effectively sam- 
ples the rest-frame UV, making these data unlikely to 
suffer biases in the estimate of host galaxy SFRs. 

The second major result of our SPS analysis was that 
galaxy mass-to-light ratios (and thus final mass esti- 
mates) were recovered without significant bias even with 
limited photometric coverage. The reason for this is re- 
lated to the aforementioned age-dust degeneracy. Both 
stellar age and foreground dust make stellar populations 
both redder and dimmer. While the ways in which the 
reddening and dimming effects are coupled (i.e. the slope 
of the respective color- luminosity relations) are not the 
same for age and dust, they are close enough that mass- 
to-light ratios are not made highly inaccurate. This 
serendipitous similarity between these effects means that 
the masses of SN la host galaxies can be fairly well re- 
covered from optical data alone. We also found that a 
single color results in significant uncertainty on the mass- 
to-light ratio of a galaxy, but this uncertainty is small 
compared to the range of masses probed by our SN la 
sample. 

We then compared our mass and SFR estimates to 
those derived using other SPS techniques. A first com- 
parison was applied to a sample of SDSS field galaxies 
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whose masses and SFRs were measured by the MPA-JHU 
group using the methods of KaufFmann et al. (2003), 
which provided the intellectual template tor our own 
methodology. We found excellent agreement with their 
values, as expected. We then used the code ZPEG, a pop- 
ular code with SN la host studies, to derive masses for 
our hosts as well as the same SDSS training sample and a 
sample of our model galaxies. We found that the masses 
showed good agreement, but the fitted SFR values from 
ZPEG converged essentially to discrete values set by the 
input galaxy evolutionary scenarios. While such a dis- 
crete set of templates is the ideal set for deriving pho- 
tometric redshifts, it provides coarse resolution for the 
detailed SFR measurements we desire. 

Furthermore, such matching to small discrete SFH li- 
braries provides an insufficient accounting for the un- 
certainty on a galaxy's SFH. While such limitations do 
not severely hinder the measurement of galaxy stellar 
masses, it provides a poor estimate of the uncertainty on 
the mass. The inability of discrete SPS models to ac- 
curate quantify the uncertainty on galaxy physical prop- 
erties was the prime motivation for our development of 
the SPS technique presented here. With the observation 
of a residual bias of SN la luminosities with host mass 
([Kelly et al]|201 0t ISul livan et al.]|2010[ [Lampeitl et al , 
2U10 " Gupta et al. 201l|), it is possible that host galaxy 



properties (especially stellar mass) may need to be in- 
cluded in SN la luminosity corrections in cosmological 
samples. In order to facilitate an appropriate accounting 
of systematic errors in this correction, we designed this 
SPS method to recover accurate measurement of the un- 
certainty on SN la host galaxy masses that accounts for 
the uncertainty in galaxy SFH. 

Finally, we comment on some of the limitations of our 
SPS methods. Our estimates of SFR and mass recovery 
efficiencies were found when fitting our model spectra 
with other models from the same library, and thus can 
be considered a minimum level of systematic error associ- 
ated with our technique. Our method is inherently tied a 
particular choice of SFHs which are meant to represent a 
reasonable prior on the likely distribution of galaxy SFHs 
in the local universe. Galaxies which are drawn from a 
different prior (i.e. if our prior is not an accurate repre- 
sentation of the true distribution of galaxy SFHs in the 
local universe) may have bia ses in the estima,tion o f their 
physical properties (see, e.g., Gallazzi & Bell 20091. Sim- 
ilarly, our prior on the reddening of these galaxies and 
our choice of reddening law may introduce some biases if 
the true distribution of galaxy reddening or the redden- 
ing law differ from the choices made here. 

While a full assessment of the possible biases result- 
ing from our choices is beyond the scope of this paper, 
we note the important fact that such biases are likely to 
affect all estimated galaxy properties in a coherent way. 
Thus a relative comparison of galaxy properties amongst 
galaxies fitted with the same SPS techniques is entirely 
valid. It is for this reason that we chose the particular 
SPS techniques presented here, because they mimic very 
closely those employed by the major galaxy studie s from 
the SPSS survey (e.g.,[K auffmann ct al. 2003; Tr emonti| 
et al.||20d4 Sahm et al. 12007) . Such self-consistency will 
also be important tor SiN la cosmological studies if SN la 
host galaxies are used for cosmological luminosity correc- 
tions. 



6.3. Comparison of SN la Hosts to Normal Galaxies 

In Sections |4] and [5] we showed that the metallicity and 
star-formation activity of SN la host galaxies is very sim- 
ilar to that of a typical field galaxy sample from SDSS. 
Not only do the mean and RMS of SN la host metallici- 
ties agree with the normal galaxy MZ relation, but indeed 
the shape of the SN la host metallicity distribution is re- 
markably similar to that of normal galaxies. An analysis 
of star formation indicators (total SFR and sSFR esti- 
mated from photometry, and Ha equivalent width and 
total flux measured from spectroscopy) in SN la host 
galaxies showed that the distribution of those quanti- 
ties in SN la hosts was similar to that calculated from 
a comparable SDSS field galaxy sample (though system- 
atic differences in SFR measurement techniques between 
the two samples - see Section |3.3| - limit the power of 
this comparison). This confirms that SNe la do not (on 
average) prefer extreme environments, but instead can 
be expected to be found in most normal galaxies. 

The most important consequence of this result is that 
studies which investigate the variation of SN la proper- 
ties with progenitor metallicity can be accomplished by 
means of measuring galaxy stellar masses from photom- 
etry, rather than the observationally expensive spectro- 
scopic measurement of galax y metallicity. In deed several 
studies ( |Neill et al.{2009; Howell et al.||2009| ) have previ- 
ously attempted to measure the decrease in S N la ^^Ni 



yield predic ted to occur at high metallicities (Timmes 
et al. |2003 ) by measuring ^^Ni yields as a function of 
host galaxy mass. Similarly, the aforementioned studies 
which discovered the SN la luminosity host bias com- 
mented on the possibility of progenitor metallicity driv- 
ing the effect. These and other works relied on the previ- 
ous untested assumption that SN la host galaxies exhibit 
typical metallicity for their given stellar mass. We have 
shown here quantitatively that such assumptions were 
(retroactively) justified, and may indeed be employed in 
the future. 

7. CONCLUSIONS 

In this paper we presented photometric and spectro- 
scopic observations of galaxies hosting SNe la discovered 
or observed by the Nearby Supernova Factory (SNfac- 
tory). Galaxy photometry in ten photometric filters 
spanning UV to NIR wavelengths were calculated for 
over 450 SN la host galaxies. Additionally, emission line 
fluxes from optical spectra were measured for over 300 
SN la host galaxies. 

We presented a Bayesian method for fltting the SN la 
host galaxy photometry with stellar population synthesis 
(SPS) techniques inspired by many of the major SDSS 
galaxy sample studies. We illustrated the vital role 
GALEX UV photometry plays in the correct estimation 
of galaxy star-formation intensity (sSFR) - especially for 
galaxies with modest to low sSFRs - in good agreement 
with t he fi ndings of previous st udies such as Neill et al. 
(12009) and [Gupta et al.| ( |201l] ). We further showed our 
methods to be consistent with other techniques for ob- 
taining galaxy stellar masses, but with improved quan- 
tification of systematic uncertainties in SPS modeling. 

Finally we used our galaxy physical parameters to com- 
pare the full sample of SN la host galaxies to a typical 
field galaxy sample from SDSS. We showed that SN la 
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host galaxies obey the fiducial galaxy mass-metallicity 
(MZ) relation with remarkable agreenient. We examined 
star-formation indicators in these SN la hosts and showed 
them to be distributed similarly as in an analogous sam- 
ple of normal field galaxies. In summary, we have demon- 
strated that SN la host galaxies are extremely normal. 
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